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r"| ' Abstract: Type I and type II seesaw contributions to the mass matrix of light neutrinos 

are inherently related if left-right symmetry is realized at high energy scales. We investigate 
implications of such a relation for the interpretation of neutrino data. We proved recently 
' that the left-right symmetric seesaw equation has eight solutions, related by a duality 

property, for the mass matrix of right-handed neutrinos M/j. In this paper the eight allowed 
^ \ structures of are reconstructed analytically and analyzed numerically in a bottom-up 

5^ \ approach. We study the dependence of right-handed neutrino masses on the mass spectrum 

of light neutrinos, mixing angle ^13, leptonic CP violation, scale of left-right symmetry 
breaking and on the hierarchy in neutrino Yukawa couplings. The structure of the seesaw 
formula in several specific 5*0(10) models is explored in the light of the duality. The 
outcome of leptogenesis may depend crucially on the choice among the allowed structures 
of Mr and on the level crossing between right-handed neutrino masses. 
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1. Introduction 



Experiments with solar, atmospheric, accelerator and reactor neutrinos have given an un- 
ambiguous evidence for neutrino oscillations and therefore for neutrino mass. This implies 
a physics beyond the standard model (SM) of particle physics, since neutrinos are strictly 
massless in that model. At the same time, neutrino mass is different from zero in almost 
every extension of the SM. Moreover, many models contain more than one source of neu- 
trino mass. A crucial issue is therefore the ability to discriminate, within a given theoretical 
framework, among these sources using neutrino data as well as other neutrino-related ex- 
perimental input. 

If lepton number is not conserved, neutrinos are Majorana particles and their mass 
originates from the dimension five operator LLtjxj) via the electroweak symmetry breaking 
(L = {vl II)^ and (f) = (f)^)'^ are the SM lepton and Higgs isodoublets, respectively). 
New physics can contribute to this operator either at tree or at loop level. We consider here 
the first possibility; the super-heavy particles exchanged at tree level can then be either 
isosinglet fcrmions Nr, that is, right-handed (RH) neutrinos [1], or isotriplet scalars 
[2]} The two cases are known under the names of type I and type II seesaw mechanism, 
respectively. Besides providing a very natural explanation of the smallness of the neutrino 
mass, the seesaw mechanism has a simple and elegant built-in mechanism for generating 
the observed baryon asymmetry of the universe - baryogenesis via leptogenesis [4] . If the 
masses of Nr and A l are much larger than the electroweak scale, then neutrino experiments 
can only probe the couplings of the effective operator. Therefore in general the data cannot 
tell type I from type II seesaw contribution to the masses of light neutrinos. Nonetheless, 
in several (possibly unified) models of fermion masses and mixing, both type I and type II 
terms are naturally present, and the knowledge of their relative size is crucial for achieving 
predictivity.^ 

In this paper we study in detail the interplay of type I and type II seesaw contributions 
to neutrino mass. We consider a theoretically well motivated connection between these two 
contributions, which stems from the proportionality of the Majorana mass terms of left- 
handed and RH neutrinos. Let the Lagrangian of the theory contain the following term: 

= -^{uICulAI + NfCNlA%) + h.c. , (1.1) 

where C = ^7270 is the charge conjugation matrix, iV£ = (Nr)'^ and A^^ ^ are neutral 
scalar fields. In the case of more than one generation, i/l and iV£ carry a flavor index and 

^We will not consider the unique remaining possibility, that is, isotriplet fermions [3], which are more 
difficult to embed in minimal extensions of the SM. 

^Even in very constrained scenarios, it is far from trivial to control the relative size of the contributions 
of the two seesaw types. This problem in the minimal SUSY 50(10) model has been extensively discussed 
in [5, 6] (for a different model see, e.g., [7]). 
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/ is a symmetric matrix of couplings. When ^ develop non-zero vacuum expectation 
values (VEVs) vl,_r, two Majorana mass terms, proportional to each other, are generated, 
one for left-handed and one for RH neutrinos. In general, a Dirac-type Yukawa coupling 
term is also present in the theory: 

CD = -yi^lCNl4P + \i.c., (1.2) 

where the neutral scalar cj)^ has a VEV v. When vji ^ v, the effective Majorana mass 
matrix for light left-handed neutrinos takes the form 

mu^vif yf^y^ ■ (1-3) 

VR 

The first (second) term on the r.h.s. of this equation is known as type II (type I) contribu- 
tion to the light neutrino mass.'^ The same coupling / enters in both contributions because 
of the assumption made in eq. (|1.1| ), which follows from a discrete left-right symmetry of 
the underlying theory. 

In fact, the seesaw formula takes the form given in eq. ( |1.3| ) in models with left-right 
(LR) symmetry above some energy scale larger than the electroweak scale [8, 9]. By this 
we mean that the gauge group contains (or coincides with) SU{2>)c x SU{2)l x SU{2)r x 
U (1) B~L and a discrete symmetry guarantees that the SU (2) l and SU (2)/j gauge couplings 
are asymptotically equal. ^ Such models incorporate naturally RH neutrinos (as well as 
isotriplet scalars) and explain the maximal parity violation of low-energy weak interactions 
as a spontaneous symmetry breaking phenomenon. Note also that LR models can be easily 
made super symmetric, and the exact R-parity at low energies can be obtained through the 
spontaneous breaking of U{1)b-l [H]- 

In this paper we perform a thorough phenomenological analysis of the seesaw formula 
( |1.3| ). A duality property of this formula, which will be at the basis of our analysis, was 
identified in a recent paper by the authors [12]. The implications of such duality in a 
specific SO{W) model are studied in [13]. 

The paper is organized as follows. In section 2 we discuss the duality property of 
the left-right symmetric seesaw. We also give the motivation for the reconstruction of 
the matrix / (and therefore of the mass matrix of heavy RH neutrinos) from the seesaw 
formula (|1.3D and study the multiplicity of the solutions. In section 3 we develop a method 
of exact analytic reconstruction of / in the cases of one, two and three lepton generations. 
In section 4 a number of numerical examples, illustrating our analytic results, are given 

^Sometimes in the literature the mechanism leading to eq. ( |l.3| ) is called type II seesaw. We prefer 
the terminology where it is called type I+II seesaw, whereas the term "type II seesaw" is reserved for the 
situation when the Higgs triplet Al is the sole source of neutrino mass. 

''classic extensions of this minimal LR symmetric setting are provided by the Pati-Salam model, based 
on SU{4)c X SU{2)l x SU{2)r = SU422 [8], and by the unified models based on 50(10) [10]. 
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and the conditions for the existence of light RH neutrinos are discussed. In section 5 we 
discuss the seesaw duahty in specific models with LR symmetry, including several 50(10) 
models. In section 6 we briefly discuss some further issues pertaining to our analysis - 
stability of the results with respect to the renormalization group evolution effects and 
baryogenesis through leptogenesis. Discussion and summary are given in section 7. The 
appendices contain further analytic treatment of relevant topics: the reconstruction of / 
in the case of antisymmetric y (Appendix A), the estimate of the relative size of type I 
and II seesaw contributions to rriu (Appendix B), and the generalization of the duality and 
the reconstruction of / in the case of a different realization of the discrete LR symmetry 
(Appendix C). 

2. Bottom-up approach to left-right symmetric seesaw 

The LR symmetric seesaw equation ( |1.3| ) relates three vacuum expectation values (f , vl and 
vr) and three n x n matrices {rriu, y and /), where n is the number of lepton generations. 
The matrix / is not directly constrained by any experiment. Moreover, being a Majorana- 
type Yukawa coupling matrix, / has no analogy with the only Yukawa coupling matrices 
that can be presently accessed experimentally (y„, y^, ye), which are all of Dirac type. 

It is therefore natural to employ the seesaw formula for reconstructing the matrix /, 
taking the quantities ruu, y, v, vl,r as input parameters. The purpose is to provide an 
insight into the underlying theory at the seesaw scale, which is deeply characterized by 
the structure of / and which is not accessible to direct experimental studies. Notice that 
eq. ( |1.3| ) is a non-linear n x n matrix equation for / and solving it in the case n > 1 is 
highly non-trivial. 

In this section we first review the experimental and theoretical constraints on the input 
parameters. Then we discuss the duality property of the LR symmetric seesaw formula 
( |1.3| ) and the multiplicity of its solutions for the matrix /. 

2.1 Parameter space 

Let us examine the possible values of the various quantities involved in eq. ( |1.3D . 

The neutrino mass matrix rrii, may be completely determined, at least in in principle, 
from low energy experiments. In terms of the mass eigenvalues rui [i = 1, 2, 3) and leptonic 
mixing matrix [/, one has 



where U depends on three mixing angles ^23 and ^13, one Dirac- type CP-violating 
phase 5 and two Majorana-type CP-violating phases p and a. From global fits of low 



my = 



U*diag{mi, m2, m^) 
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energy experiments [14, 15] one finds 

Amli = ml-m\ = {J.dtH) x 10"^ eV^ , Am^^ = \ml-ml\ = (2.6±0.6) x lO^^ eV^ , 

(2.2) 

ei2 = {S3+l)\ 023 = (45 ± 10)° , 013 < 12°, (2.3) 

where the best fit values and 3a intervals are taken from the last update of [14] (v5). 
Cosmology yields the most stringent upper limit on the sum of the neutrino masses, mi + 
1712 + TTT'S ^ 0.4(0.7) eV at 99.9% C.L. with(without) including Lyman-a data in the fit 
[16]. The values of the CP-violating phases are completely unknown at present. 

As far as the vacuum expectation values are concerned, we choose the convention of real 
and positive v, vl and vr (this can always be achieved by redefining the phases of ul, 



and y in eqs. (1.1) and ( |1.2| )). The electroweak symmetry breaking parameter t; ~ 174 GeV 
is accurately known, while only weak constraints are available for the triplet VEVs 
No lower bound on f^. exists, and actually conditions are known under which Al does not 
acquire any induced VEV. In contrast to this, vl is bounded from above by its contribution 
to the p-parameter (Ap ~ — 2w£/f^), so that precision electroweak measurements imply 
Vl ^ GeV) [17]. Since we do not know the scale for the onset of LR symmetry, there is 
essentially no upper bound on vji; we will be assuming in the following vr < Mp\. The 
value of vr is bounded from below by the non-observation of RH weak currents (vr > TeV) 
[17]. It is useful to translate these bounds into bounds on the parameter x = vlvr/v^ . If all 
dimensionless couplings in the scalar potential are of order one, its natural value is x ~ 1. 
However, the scalar sector may well be very complicated and the potential may depend on 
various mass scales which are not constrained a priori. Therefore, x should be considered a 
free parameter in a model independent analysis. Using the above experimental constraints, 
one finds < a; < 10^^. 

The Dirac-type neutrino Yukawa coupling matrix y is not directly measurable, at least 
at present. Nonetheless, in LR-symmetric models and in their embedding in models with 
partial or grand unification, y is usually related to known quark and/or charged lepton 
Yukawa couplings. A detailed discussion of the structure of y in some of these models will 
be given in section ^. Here we constrain ourselves to a few examples. In the minimal LR 
symmetric model, in the supersymmetric case one has y = tan f3ye, where j/e is the Yukawa 
coupling matrix of charged leptons and tan [3 is the usual ratio of VEVs. In the minimal 
Pati-Salam model, one obtains instead y = yu, where is the Yukawa coupling matrix of 
the up-type quarks. This relation also holds in the 50(10) model where the only Higgs 
multiplet contributing to the fermion masses is in the fundamental representation 10 h- 
More complicated (and realistic) examples will be discussed in section ^. 

One should bear in mind, however, that the nature of the high energy theory is es- 
sentially unknown and therefore the Yukawa coupling matrix of neutrinos may in fact be 
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very different from those of charged fermions. In particular, the neutrino flavor sector 
may have radically different symmetries. As an interesting example, one may consider 
y = diag{a,b,b) [18], which is motivated by the maximal or nearly maximal 2-3 mixing 
observed in the atmospheric neutrino oscillations. It should be also noted that, in the 
presence of low energy super symmetry, lepton flavor violating processes can be directly 
sensitive to the neutrino Yukawa couplings. For example, in constrained mSUGRA seesaw 
models one may need yij <C 1 for some i,j, in order to suppress processes like t ^ fi^ [19]. 

This brief survey illustrates that although certain choices of y may be well motivated, 
this matrix is in fact model dependent. In the following we will keep the form of y as 
general as possible, so as to maintain the generality of our results within the chosen LR 
symmetric framework. 

2.2 Seesaw duality 

Before discussing the general case, it will be useful to find the matrix / approximately in 
the limits when either of the two contributions to rrii, in eq. (1.3) dominates (these limits 



are well defined as long as all the matrix elements in rrii, are dominated by the same type 
of seesaw). In the case of dominant type I seesaw, one obtains 



(2.4) 



m = , X = — 2~ • l^-^j 



// = y'^m ^y -\ — 7^{ym ^y)'^m ^(ym ^y) 

X x^ 

Am'^^ymT^y)^ (v^ m~^y + ym~^y^)(m^^ym^^y) + . . . 

x-^ 

where we used the notation 

nil, vi 

m = q' = 

Vl W 

In the case of dominant type II seesaw, one finds 

fii = m + -ym^^y'^ \{ym~^y)nr^{ym^^y)'^ + ... . (2.6) 

X 

In this latter case my ~ vl/u, so that low energy neutrino data allow one to reconstruct 
directly the mass matrix of RH neutrinos M/j ~ vrJij. In particular, the spectrum of the 
heavy neutrinos coincides (up to an overall factor) with that of light neutrinos [20]. 

In general, there is no a priori reason to expect that one type of seesaw dominates. 
If both contributions to rrii, are comparable, the solution of eq. ( |1.3D cannot be obtained 
expanding around the purely type I or II solution and a different approach is needed. 

The LR symmetric type I+II seesaw formula in eq. ( [1.3| ) has the following intriguing 
duality property. Suppose that a matrix / solves this equation. Then one can verify that 

f = m-f = --yf-V (2.7) 

X 



-6- 



is also a solution, provided that the matrix y is invertible and symmetric (or antisymmetric). 
These conditions on y turn out to be both necessary and sufficient for / to solve eq. (|1.3|). 
In particular, they guarantee that the matrix / is invertible, as it must be in order to 
satisfy eq. (|1.3|). It is easy to see that the duality operation is closed, i.e. the dual of / 
coincides with /. 

An important example of duality is provided by the solutions for / in the case of one 
seesaw type dominance, given in eqs. (|2.4|) and (2.6): if y = ±2/"^, one immediately sees 



that // and /// are dual to each other. This allows us to establish an important result: if 
there is a solution of the seesaw equation with dominant type I seesaw, then there is also 
an alternative one, with dominant type II seesaw, and vice versa. More generally, we will 
show that the LR symmetric seesaw equation always has a pair of dual solutions which 
reduce to // and /// when the input parameters satisfy certain conditions. 

In the following we will focus on the case where y is symmetric. This is true in 
models with the discrete LR symmetry ul Nf^ (see eq. (|1.2| )), including the minimal LR 



symmetric model and 50(10) models where only the couplings to 10// and 126// Higgs 
multiplets contribute to y (for further details see section ^). The case of antisymmetric y 
will be discussed in Appendix |A[ When y is symmetric, one can write y = UydU^, where 
yd is diagonal, real and positive and C/ is a unitary matrix. (We consider nii, in the basis 
where the mass matrix of charged leptons is diagonal, so that U describes the mismatch 
between the left-handed rotations that diagonalize the Yukawa couplings of charged leptons 
and neutrinos.) As a consequence, eq. ( |1.3| ) can be rewritten as 

U^m.U* = VLiU^fU*) - ^ya{U^ fU*)-^yd . (2.8) 

VR 

Therefore, one can always work (and we will) in the basis where y is diagonal, by redefining 
{Wm,yU*) —>■ niy and {W fU*) — > /. If f7 ~ 1, the input from the low-energy neutrino 
data will still (approximately) determine the left-hand side of eq. ( |2.8| ). Note that U is the 
leptonic analogue of the CKM mixing matrix in the quark sector, where mixing is known to 
be small. Therefore the condition ~ 1 may be motivated by a quark-lepton symmetry. 

2.3 Multiplicity of solutions 



We have found that, if a matrix / solves the seesaw equation (L3), so does / = m — /, i.e. 
this equation does not have a unique solution. This is hardly surprising, as eq. (|1.3| ) is a 
non- linear matrix equation for /, or, equivalently, a system of non- linear coupled equations 
for its elements fij. From the duality property of eq. ( |1.3D it immediately follows that the 
number of its solutions must be even. We shall now show that for n lepton generations the 
multiplicity of solutions is 2". ^ 



^This result was first obtained in [12], tiiougli in a less straightforward way. An alternative derivation 
can be found in [131. 
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Let us introduce the matrices m and / through the relations 



1 



,1/2 ^,,1/2 



m = —^y ' my 



f = ^y'"h''' 



(2. 



where y^^'^ satisfies (y^^^)^ = y- From eq. (1.3) with symmetric and invertible y one then 
obtains 

m = /-/-!. (2.10) 

The duahty we discussed above is especially clearly seen in this equation, where it corre- 
sponds to the invariance with respect to / ^ —f~^ {= m — f). Multiplying eq. ( p.lO| ) by 
/ on the left or on the right, we find that / satisfies 



f-mf -1 = 



(2.11) 



and that [rh, f] = 0. The commutativity of fh and / allows one to find a formal solution 
of eq. ( 2.11| ). Indeed, let us write / in the form / = m/2 + R, where the matrix R is yet 



to be determined. Substituting this into eq. ( 2.11 ) in which the term rhf is rewritten as 



(l/2){m, /}, one finds that the matrix R satisfies i?^ = ^ + 1, so that finally one obtains 



~ m / 



(2.12) 



Since any non-singular n x n matrix with non-degenerate eigenvalues has 2" square roots 
[21], eq. (|2TC| ) (and so eq. ( |l.3D ) has 2" solutions. Obviously, if Rq is a square root of R^, 
so is —Ro', therefore the solutions in eq. ( ^l2|) form 2"-i dual pairs with f + f = fh. 



3. Analytic reconstruction of the coupling matrix / 

In this section we will solve analytically the seesaw equation for / in the cases of one, 
two and three lepton generations. For one generation, the analysis is straightforward; 
for more than one lepton flavor, the effects of mixing complicate the reconstruction of / 
considerably and, for this purpose, we shall need to develop new algebraic techniques. We 
will provide general solutions for the case of symmetric y, where the seesaw duality occurs 
(the generalization to antisymmetric y is given in Appendix 

We will also identify the criteria to quantify the dominance of one or the other type 
of seesaw. In the presence of mixing, this identification will turn out to be a subtle issue, 
whose details are given in Appendix ^. 

A short account of the main results of this section was published in [12]. An alter- 
native approach for the analytic reconstruction of / was developed in [13], for the case of 
symmetric y. 
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3.1 The case of one lepton generation 



In this case my, y and / are merely complex numbers. Eq. (O), being quadratic in /, has 
two solutions, which we denote f±: 



1 



A 2 2 \ 1/2' 

VR 



l±{l + d) 



1/2 



where 



d 



2, ,2 



(3.1) 



(3.2) 



Obviously, the solutions f± are dual to each other: m = /+ + /_. 

At this point it is useful to state our convention for the assignment of the complex 
phases. The freedom to rephase the fields in eqs. (LI) and ( |L^ ) allows us to take vl, 



vr as well as rriy real and positive. Next, we define = |ype*-*^, so that argd = x ^iid the 
phase of / is then determined by eq. (|3.lD . 

Type I and II contributions to rrij^ are defined as 



m„ 



vr! 



m 



II 



(3.3) 



Their relative size is determined by the value of d: 



ii/i 
r_|_ 



m 



II 

u± 



[i±{i + dY/'^Y 

d 



m 



(3.4) 



We identify three possible physical regimes: 



(i) \d\ « 1 

(ii) \d\ 1 

(iii) > 1 



1 772"^ I ~ 1 771"^"^ I 



-m 



II 



and \mlL\ <C \ml_ 



single type dominance 
hybrid seesaw 
cancellation regime 



Let us discuss these cases in turn. 

(i) The dominant seesaw type is I (II) in the case of the /_ (/+) solution. For \d\ <C 1, 
eq. ( p. ID becomes 

9 9 9 9 

J- ~ , /+ ~ \ • (3-5) 

VL VRin, 



VRTTly 



Therefore a value |d| ^ 1 implies that one type of seesaw is dominant, but it does not 
determine which one. Thus, both type I and type II dominance limits correspond to the 
same condition \d\ <^ 1. This may look counter-intuitive, as these two limits correspond to 
apparently opposite conditions \ml\ » and \ml\ <C \ml^\. However, it is easy to see 

that, when expressed in terms of the "input" parameters only, i.e. when the corresponding 
solutions for / are substituted, both conditions reduce to |d| <C 1. 
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eV 

1 r 




d 



Figure 1: Moduli of the solutions VLf± of the seesaw formula, as functions of We chose 
■nil, — V ^'rriatm — 0-05 eV. The solid, dashed and dotted curves correspond to argd = 0, 7r/2 and 



TT, respectively. Recall that rrij/^ — t'L/± and 



= — m 



A solution / is meaningful (perturbative) only if |/| < 1: for ij^^ <C rriu, the solution /_|_ 
violates perturbative unitarity and should be discarded, so that only /_ is viable (dominant 
type I). A necessary condition for both solutions to be perturbative is vr > v^\y\^ /my, so 
that either type I or type II contribution to rrii, can dominate. Notice that for jdj ^ 1 
the expansions ( |2.4| ) and (^^) apply for /_ and /+, respectively; eq. ( |3.5| ) just gives a 
simplified version of such expansions. 

(ii) For ~ 1 both seesaw types give sizable contributions to rriy. Notice that when 



d = -I (i.e., AvLv'^\y\'^ 



VRrriy and x = argy^ = vr), the equality ml 



m 



II 



my/2 is 



realized. This degeneracy point corresponds to the absence of CP violation. ^ 

(iii) The case |(i| ^ 1 corresponds to a cancellation between the two seesaw contribu- 
tions to my. To bar a too strong ("unnatural") cancellation, one may demand \my^^\ < 1 
eV. Eq. ( |3.3| ) then implies a lower bound on the scale vr of SU {2)r symmetry breaking: 
VR > v'^\y\'^/i\f \ -l eV) ^ |yV/|3 • 10^^ GeV. If one makes the additional ("natural") 
assumption x = vrvl/v'^ ~ 1, eq. ( |3.3D will also imply vr > u^|/|/eVss |/|3 • 10-*^^ GeV. 

For illustration, let us assume that the light neutrino mass is given by the atmospheric 
scale, my = 0.05 eV. The moduli wl|/±| of the solutions of eq. (3.1) as functions of \d\ are 
then shown in fig. ||. The three regimes (i), (ii) and (iii) correspond to the left, central 
and right regions of the plot, respectively. The value of |/+| (|/_|) is a slightly decreasing 
(increasing) function of arg d = x between and tt, as shown in the figure. For X = ^ 
|ci| > 1 the moduli of the two solutions coincide. 



®Note that in general CP violation is present in the one-generation case due to the presence of the Higgs 
triplets. Although it does not manifest itself at low energies, it can lead to a successful leptogenesis (more 



on this in section 3.2) 
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GeV 




Figure 2: The mass Mii± — w_r|/±| of the right-handed neutrino for the two solutions f± of the 
seesaw formula. We chose to„ = -y/Arof^ = 0.05 eV and \y\ = 1. For fixed Mr scales as 
The solid, dashed and dotted curves correspond to argc? — 0, 7r/2 and tt, respectively. 



The mass of the RH neutrino is given by 

2y2[l±(l+d)i/2] 



MR±=VR\f±\ 



d 



(3.6) 



In fig. |2| we plot Mr as a function of \d\ for = 0.05 eV and \y\ = 1. The mass of 
the scalar triplet Ma is not directly related to the parameters in the seesaw formula. One 
expects Ma ~ /vl ~ vr, but very different values are also naturally possible, depending 
on the details of the mechanism which induces v^. 

Consider now the case of n lepton generations with no mixing (which means that / 
and y are diagonal in the same basis). This case can be described as n replicas of the 
one-generation case, each with two solutions {i = l,...,n). For each generation, the 
relative size of type I and II contribution to the light neutrino mass rriui depends on the 
dominance parameter 

= — J' , i = l,...,n. (3.7) 

For a given value of vi/vr, one can have \di\ <C 1 (single type dominance) for an i-th 
generation, while at the same time the hybrid seesaw or cancellation regimes may be 
realized for j ^ i- We will see that in the presence of flavor mixing it is much less trivial 
to assess the relative size of the two seesaw types. 

3.2 The case of two lepton generations 

In this case in the flavor basis, where the mass matrix of charged leptons is diagonal (say. 
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mi = diag{m^,mr)), the seesaw formula ( |1.3| ) can be explicitly written as 

^ / /22 /23 \ ]_ / 2/^2 \ / /33 -/23 \ / y,i2 yr2 \ 

V/23 /33 / \ yT2 yr3 J \-f23 /22 / 1 2/^*3 yr3 / 




where F = det f = 722/33 — /23- The matrix equation (|3.8|) is equivalent to the system of 
three coupled non-linear equations for /22, /33 and 723: 

2;-^(/22 - rUf^f,) = /33y^2 - 2/23^^22/^3 + /22y^3 ' 

2;-^(/23 - m^r) = /33yM2yT2 " /23(y/i3yr2 + yti2yr3) + /22y/i32/T3 , (3-9) 
a;-^(/33 - rrirr) = f33yr2 " 2/23yr2yr3 + /22y?3 • 

To solve the system (|3.9|) , we use the following procedure. Let us define /' = f/^/X, 
m' = m/\/A and y' = y/\/A, where A is an as yet arbitrary complex number. The scaling 
law was chosen in such a way that in terms of the primed variables the system of equations 
for flj has the same form as eq. ( |3.9D . Next, we fix the value of A by requiring F' = 
detf = 1. The system of equations for f-j then becomes linear and can be readily solved. 
Expressing the primed variables back through the unprimed ones and substituting them 
into the condition F'(X) = 1, one obtains a 4th order polynomial equation for A. In general, 
it has four complex solutions Aj (i = 1, . . . , 4), leading to four allowed matrix structures /j. 
Notice that Fj = det /i = Aj. This procedure proves that eq. ( |3.8[ ) has four solutions for a 
generic structure of y and, therefore, it generalizes (for the case n = 2) the proof given in 



section 2.3 for symmetric y. 

We now present the general analytic solution in the case where the matrix y is invertible 
and symmetric, so that the duality holds. As follows from eq. ( |2.8D , when y = y'^ we can 
choose the basis where y is diagonal: 2/^3 = yr2 = 0, y^2 = 2/2, 2/t3 = 2/3- With a little 
abuse of notation, we will still denote the matrix elements of m and / in the new basis as 
rriap and fij. The system (3^) can be linearized as described above and is easily solved: 

^ _ xX I xXm^^ + yjmrr m^AxX- 



xXy - y^yi \ ... xXm^r + yl 




m 



where A is a solution the following quartic equation: 

[{xXf - 2/22/3] -X [det m{xX - y2y3fxX + {m^^ys + mrTy2f{xX)'^] = . (3.11) 

Taking the determinant of the equality f = m — f = —yf^^y/x, one obtains 

x^AA = x^F-F = yjyl , (3.12) 

where F = det /. With the help of this relation, it is straightforward to check that the four 
solutions fi defined by eqs. (|3.10| ) and (|3.11| ) form two dual pairs. Actually, the duality 
makes it easy to express the four solutions of eq. ( p. 11 ) in the closed form as follows: 

xXi = - X det m + r± ± 2{x det m)^ + 4kx + 2r±x det m , (3.13) 
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where 



m 



ys + 2m^^y22/3 + rn^ry2 > r± = ±J{x det mf + ikx + 16y|y| . (3.14) 



One pair of dual solutions corresponds to and the other one to r_; within each pair, 
a solution is distinguished from its dual by the sign in front of the radical in eq. ( 3.13| ). 



Summarizing, eqs. ( |3.1[1| ) and ( 3.13| ) define explicitly the four solutions fi as functions of 
the input parameters m = m^/vL, x = vlvr/v^ and y. 

We can now study the structure of fi in any regions of parameters of physical interest. 
To start with, let us consider the limit 2/2 0. This can be justified if neutrino Yukawa 
couplings are related to those of charged fermions, so that 2/2 ^ ^ 1- Then eq. ( 3.13| ) 
can be expanded as 

xXifi = ^ \ + c(y|) , = xXi^ = -^-^ , (3.15) 

where = ziz^(x det m)^ + Axy'^vn?^^. The solutions for / take the form 

/i,3 = m- /2,4 , /2,4 = y2^^ , I + 0{yl) . (3.16) 



V " a:Ai,3 

Since A2,4 = det /2,4 are proportional to in the case of the solutions /2,4 one RH neutrino 
mass becomes much smaller than vr for very small 7/2- Notice that, even though A2,4 go 
to zero for y2 0, the matrices /2,4 are finite and invertible (and therefore acceptable 
solutions) for any 7/2 7^ 0. 

Let us define, in analogy with eq. (|3.2|), a dominance parameter that controls the 
relative size of type I and type II seesaw contributions: 



VL Av'^yl _ 4:y? 
vr {my)1^ xm? 



d^^Zl^ = ^. (3.17) 



Obviously, other analogous parameters can be defined by replacing m-j-r and/or y^ with the 
other entries of m and y. A detailed discussion of this issue is postponed to Appendix 
here we just study the main features of the dependence of the solutions on d. Fig. |^ shows 
the dependence of the values of |xAi| on \ d\ for a specific realistic set of the input parameters. 
In the region \d\ <^ 1, there is a solution |xAi| ^ \yiyj\, which leads to /i ~ m (dominant 
type II seesaw), while for the dual solution, one has |xA2| <C \yiyj\ (dominant type I seesaw). 
In general, the other pair of dual solutions corresponds to hybrid seesaw. Only if the value 
of det m is strongly suppressed, one finds fs ~ /i (dominant type II seesaw) and /4 ~ /2 
(dominant type I seesaw). In the region \d\ ~ 1, all four solutions are of hybrid type. 
Finally, for \d\ ^ 1 a cancellation between type I and II seesaw contributions to rrijy occurs. 
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|XA2| 



Figure 3: The four solutions xXi of eq. ( 3.13 ) as functions oi \d\ = 4y|/(a;TO^^). We chose 



"1^^ — TTLrT = (0.02 eV)/vL and m^r = (0.03 eV)/vL, so that 623 = 7r/4 and Am23 = 2.4 • 10 
eV. We also took yz — I and 2/2 = 10^^. 



GeV 



GeV 





0.01 0.1 1 



0.01 0.1 1 



100 1000 lOOOi 



Jd| 



Figure 4: The masses M2,3 of the RH neutrinos as functions of \d\ for the same set of input 
parameters as in fig. ^ Thick solid (dashed) curves in the left panel correspond to the solution /i 
(/2), thin solid (dashed) curves in the right panel, to the solution /s {/a). 

In fig. ^ we plot the values M2,3 of the two RH neutrino masses (that is, the eigenvalues 
of Mr = vrJ) versus \d\ for the same set of the input parameters as in fig. ^. In order 
to assess the size of M2,3, it is useful to define a "seesaw scale" Mg = /rri'^, which for 
m3 = 0.05 eV is equal to 6 • 10^*^ GeV. Since for \d\ <^ 1 the solution fi (/2) corresponds 
to type II (I) seesaw dominance, one has M2,3 » Mg for /i and M2,3 ~ y^s^s /s) a-s 
shown in the left panel of fig. Q For the solutions f^^4, the situation is intermediate and 
one finds M3 3> Mg, M2 ~ y2 3Mg, as shown in the right panel of fig. |^. For \d\ » 1, the 
asymptotic values of M2 are the same for all four solutions; the same is also true for M3, 
with M2/M3 ~ 2/2/2/3- For the solution /i, one observes a level-crossing of the two RH 
neutrino masses for \d\ ~ 20. One could further characterize the sector of RH neutrinos by 
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studying the dependence of their mixing on d. 



3.3 The case of three lepton generations 

In this case, the symmetric matrix / contains 6 independent elements fij. The LR sym- 



metric seesaw relation (y_^) can be written as a system of 6 nonlinear coupled equations for 
these elements, analogous to that in eq. ( |3.9D . Here we solve this system analytically for 
the case of a symmetric neutrino Yukawa coupling matrix y, for which the duality holds. 
Using the freedom to work in the basis where y is diagonal, y = diag{yi,y2,y3), one can 
write this system as 

xF{fij - rriij) = yiyjFij , (3.18) 

where F = det f and 

_ 1 

-^ij = '^^ikl^jmnfkmfln ■ (3.19) 

The six corresponding equations for the dual matrix f = m — f are 

xF{fij - rriij) = -xFfij = yiyjFij , (3.20) 

where F = det / and 



Fij — ^ikl^jmnfkmfln — F^j Tij -\- Mij , 
_ ^ _ 1 

Fij = (^ikl^jmnfkm'm-ln > -^ij ~ '^^ikl^jmn^km^ln 



(3.21) 



Taking the determinant of the equality / = —yf ^y/x, one proves immediately that 

x^FF = -yhlyl , (3.22) 

which is an analogue of eq. ( 3.12| ). 



To solve the system ( ^.18 ) for /, we make use of a procedure similar to that employed 
in the 2-generation case. Let us define /' = //A^/^, m' = m/\^^'^ and y' = y/X^^^ where 
A is determined from the condition F'{X) = det /'(A) = 1. At this point the left-hand 



(right-hand) side of equations (3.18), written in terms of the primed variables, becomes 



linear (quadratic) in f^^. Using the duality (eqs. (|3.20|) - (|3.22| )), the quadratic terms, F/^, 



can be rewritten as 



4 = ^#^4;-^+^- (3-23) 



This allows one to linearize the system for 



- {y'lVWsfUl, - x^rn[^ = x^yWjiTl, - Mi^) . (3.24) 

After some algebra, the solution for / = X^^^f can be written in a rather compact form 
(to be compared with eq. ( 3.1[1| )): 



_ A2 [(A2 - y2)2 _ det m + Y^S] rrnj + A (A^ - F^) - Y^X'^iX'^ + Y^)Sij 
•''^ ~ (A2 - y2)3 _ y2^2(^2 _ y2)5 _ 2y2^3 det m ' 

(3.25) 
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where 



VkVl 



Ai 



E 

k,l=l 



The value of A is determined from the equation 



F(A)^det/(A) = A, 



mkix\ 

ykyi J 

(3.26) 



(3.27) 



where /(A) was defined in eq. ( ^.25 ). Eq. ( 3.27 ) turns out to be an 8th order polynomial 
equation for A, ^ which has in general eight complex solutions. Defining 



/ ykyiM, 



kl 



— ' V X 

1=1 ^ 



(3.28) 



this equation reads 



[(A2 _ y2)2 _ Y^X^sf - A2(A2 + Y^fA - Y^X^detn 
-A [A^ + ^2^2(^2 _ y2) (5 + 5) _ y6j det m = , 



(3.29) 



to be compared with eq. (3.11). The eight solutions of eq. ( 3.29] ) form four dual pairs 
Aj, Aj [i = 1, . . . ,4). The duality is described by the relation AjAj = —Y"^, which follows 
from eq. ( |3.22| ). As a consequence, eq. ( |3.29| ) can be rewritten as 



= n(A-A,)(A+n=n(A' 

i=l ^ i=l 



y2 

ZiX - , z, = Xi- — . (3.30) 

A,; 



By comparing eqs. ( p. 29 ) and ( p. 301 ), it is easy to verify that Zi are the roots of the following 
quartic equation: 



- det mz-^- i2Y^S + A)z^- Y^ (8 + S)detmz + Y' [Y^S' -4A- (det my] = . (3.31) 

Thus, although the general order 8 algebraic equation does not admit analytic solutions, 
the equation for A can be solved in radicals since, due to duality, it reduces to a quartic 
equation in z. For the latter a general solution in radicals is known, though complicated. 
The 8 solutions for / are obtained by plugging into eq. (3.25) the values of A given by 



A.(A,) 



Zi ± Jzf + 4y2 



1,...,4), 



(3.32) 



where the sign in front of the radical distinguishes between Aj and Aj 



To recognize that, one has to use duahty relations to identify and cancel a common polynomial factor 
in the numerator and denominator of F{\). 
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Given the complicated algebraic form of the solutions for /, it may seem a hopeless 
task to quantify the dominance of one or the other seesaw type in the eight cases. However, 
it turns out that the techniques developed for analyzing this issue in the one- and two- 
generation cases can be generalized here. The full details of the analysis will be given in 
Appendix |^. Here we only describe the generic case, in which the relative size of type 
I and type H contributions does not change much from one element of the matrix m to 
another and, in addition, no special cancellations among the entries of m, such as leading 
to |detm| <C |mQ,^p, occur. Then, in terms of the dominance parameter \d\ defined in 
eq. ( p.lTp , the classification goes as follows: for |(i| ^ 1 there is one pair of dual solution 
with one seesaw type dominance, while the other three pairs correspond to hybrid seesaw. 
All three RH neutrino masses are generically larger than /rrii in the solution with type 
n dominance, only two of them are larger than /rrii in three of the hybrid solutions, 
only one mass satisfies this condition in the corresponding three duals, and finally all RH 
neutrino masses are of the order of /rrii in the case of dominant type I seesaw. When 
\d\ > 1, all 8 solutions for / lead to hybrid seesaw. 

In the next two subsections we will specialize our general analytic solution to some 
physically interesting limits. The reader more interested in numerical examples may pro- 
ceed directly to section 



3.3.1 Limit of hierarchical neutrino Yukawa couplings 

A considerable simplification of the general solution for / occurs when type I contributions 
to the elements of rriy, proportional to yi, are negligible, i.e. in the limit yi ^ 1. In fact, 
this limit is physically well motivated, in view of the tininess of the Yukawa couplings of 
the charged fermions of first generation. Strictly speaking, for yi = the matrix y is not 
invertible, so that the duality of solutions for / does not hold. Nonetheless, eq. ( p.25| ) is 
valid for any small but non-zero value of yi and therefore can be expanded in powers of 
yi. Assuming that A is finite in the limit yi — > 0, one finds 



nip 




me. 



\ 



2/3 



nippm 



m 



xX 



11 1 

2/2 2/3 "^ee 



1 + 



yiysrrie 



xX 



xX 



111 

yiyi^ee 



\ 



{xXf 



+ 0{yi) 



J 



(3.33) 

The e-row of the mass matrix of light neutrinos m directly determines the first row of / 
(pure type II seesaw), whereas the 2 — 3 sector of / depends on the interplay of the elements 
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of the e-row and /ir-block of m, as well as on the values of 1/2, ys and xX. 

Consider now the behavior of A in the limit yi — > 0. For yi = 0, the 8th order 



polynomial equation (|3.29[) reduces to 



{ [{xXf - ml^ylylf - x [det m(xA - meeyiyzfxX + (M222/2 + M^^^y^if {xXf] } = . 

(3.34) 



The zeros of the term in curly brackets in eq. (3.34) are (to leading order) the four 
solutions A which are finite for yi — > 0. The four duals A = — yfylyi/l^^'^) vanish in this 
limit, so that for them eq. ( |3.33| ) does not hold. However, the corresponding solutions can 
be obtained by duality: f = m — f with / given in eq. ( |3.33| ). By expanding directly 



eq. ( 3.25| ), it is easy to check that /n ~ yf, /i2,i3 ~ yi, while the other entries are finite. 



As a consequence, one RH neutrino mass is much smaller than vr since it is proportional 
to yf. 

When yi = 0, one has only four (instead of eight) solutions, since their duals become 
singular. From the comparison of the term in curly brackets in eq. ( 3.34| ) with eq. (3.11), 



a strong analogy with the pure two- generation case becomes evident. In fact, a different 
duality among the four remaining solutions is present: if A / satisfies eq. ( p.34 ), also A = 



yiyi^^ee/l^^-^) does, and it corresponds to / = m — /, where rha/s = "map + niea^ep/n^ee- 
There are two pairs of such solutions. 

A very simple yet non-trivial scenario corresponds to the case when type I contributions 
to nil, proportional to 7/2 are also negligible (7/2 ^0). In this limit the relation fap = rriap 
(pure type II seesaw) holds for all the entries of / but /sa [20, 22]. There are two solutions 
for /33, defined by 

(/33 - m^rf + ^(/33 - rurr) - ^ = . (3.35) 

M33 X 

In particular, this equation admits I/33I l^n-rrl ~ l^^/i^tl) that is, / with hierarchical 
structure can lead to rriy with large 2 — 3 mixing. 

3.3.2 Limit of hierarchical Hght neutrino masses 

When the mass spectrum of light neutrinos has normal (inverted) ordering, the lightest 
neutrino mass mi (ms) can be negligibly small. In this case the determinant of m = rriu/vL 
vanishes. For detm = (more generically, for | detm| <C \Y\), eq. ( p. 29 ) reduces to a pair 
of quartic equations: 

(A^ - y2)2 _ [±^a(A2 + y2) + y2^A2] = . (3.36) 

If Xi is a solution of the "sign +" equation, then also Y'^/Xi is, while — Aj and —Y'^/Xi are 
solutions of the "sign — " equation. Therefore, two pairs of dual solutions for A are equal 



in absolute value and opposite in sign to the other two pairs. Moreover, eq. (3.36) has the 
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same form as eq. ( ^.11 ), if one makes the identifications 

y2 ^ ^213 ^ det m , S ^ 



xk 

2 2 ' 



(3.37) 



where the quantities on the right hand sides refer to the two-generation case. Therefore, 



the quartic equations ( |3.36D can be exphcitly solved and analyzed as in section p^ . 

Since the solar mass squared difference is much smaller than the atmospheric one, in 
the case of the normal hierarchy one can neglect, in first approximation, both mi and m2- 
In this limit becomes a rank-1 matrix, and one can write 

/ ab ac\ 



m 



ab 6^ be 
ac be 



(3.38) 



Almost maximal atmospheric mixing implies 6 ~ c as well as ^13 ~ a/b. The solar mass 
scale is zero in the limit of eq. ( p.38| ), so that the 1-2 mixing is undefined. Notice that 
plugging eq. ( 3.38| ) in eq. (|3.26|) , one finds Aij = and Sij = rriijS. As a consequence, the 
expression for / given in eq. ( p. 25 ) becomes singular in this limit, since / oc m and thus 
is not invertible. One therefore has to resort to a different approach in order to find /. It 
turns out that in this limit the seesaw equation ( p. 18 ) has an infinite number of solutions. 
To illustrate this phenomenon, consider for simplicity the case yi = 0. Eqs. (|]3|) and 
(3.34) cannot be used since they would yield f = m, an unacceptable result when the 
matrix m is not invertible. A direct solution of eq. ( p.l8| ) gives in this case 

^a^ ab ac ^ 

ab b'^ + y2Cosa be + ^2/2^3 sin a 



f 



^ac bc + ^2/22/3 sin a c^-ys cos a j 



where a is an arbitrary complex number different from zero. (For a 
ambiguity appears: 



(3.39) 



a different 



/ 



/ ab ac 
ab 6^ ± 2/2 be 
be -2 



\ 



(3.40) 



Y ac be ± 2/3 / 

where the two "it" signs are uncorrelated so that, for this value of q, there are four 
solutions.) Summarizing, when the matrix m is of rank 1 there is an infinite number of 
solutions for /. 

4. Numerical examples 

We have shown that, for arbitrary values of 2/1,2,3; vl,r and of the elements of m,^, one 
generically finds eight solutions for the matrix /, defined by eqs. ( 3. 25] ), ( |3.32| ) and ( 3.31| ). 
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The allowed values of the input parameters were discussed in section 2A. Here we present 
several numerical examples which are realistic, in the sense that they reproduce the ob- 
served neutrino oscillation parameters and respect all the other experimental constraints. 
This will enable us to identify the most interesting features of the allowed structures of /. 

1. As a first realistic numerical example that we will use as a benchmark point in the 
space of the input parameters, let us take 



m 



VL 



( 0.1 -0.1 
0.1 0.55 0.45 
y-O.l 0.45 0.55 



(4.1) 



which corresponds to the tri-bi-maximal mixing (tan^ 023 = 1, tan^ Q\2 = 1/2, tan^ ^13 = 0), 
no CP violation and Am^^^/Am^j^ ~ 0.031, in agreement with all the current data. 
The eigenvalues of m are (—0.1, 0.2, 1), so that the spectrum has the normal hierarchy 
and the lightest neutrino has CP parity opposite to that of the other two. The above 
choice of m fixes also the value ~ Am^^^ 0.05 eV (see eq. ( p. 21) 



^ ^ We also take 

X = vivji/v^ = 1, which is a natural value if the dimensionless parameters in the scalar 
potential are all of order one. This determines vr Q ■ 10^^ GeV. Finally, we take the 
hierarchy among the eigenvalues of y to be slightly weaker than that for Yukawa couplings 
of the charged fermions: yi = 10"^, 1/2 = 10~^, 7/3 = 1. We neglect possible (CKM-like) 
rotations between the flavor basis and the basis where y is diagonal. 
With these choices, the solutions of the LR seesaw formula are 

/ -0.001 0.10 -0.14 \ / 0.001 -0.005 0.04 \ 

A ~ ... -0.006 -0.04 



-0.001 0.10 -0.14 
... 0.56 0.49 
0.88 



/2 



/3 



-0.01 0.11 -0.04 
... 0.55 0.44 

-0.88 

0.02 0.07 -0.02^ 
... 0.61 0.30 

1.58 

0.01 0.08 0.08 
... 0.60 0.25 
-0.19 J 



f2 



h 



\ 



-0.33 

006 -0.008 -0.06 
.. -0.004 0.01 
1.44 

-0.02 0.03 -0.08^ 
... -0.06 0.15 
-1.03 

-0.01 0.02 -0.18 
... -0.05 0.20 
0.74 



(4.2) 



The rounding off in the numerical values of fij is chosen so as to clearly illustrate the 
matrix structure of the solutions. We have checked that for the high precision solutions 
the duality relation fi + fi = mis satisfied very accurately for each dual pair. The structure 



( |4.1| ) of m is recovered by plugging fi or fi back into eq. (|1.3| ) , but very accurate values of 
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their entries (not shown in eq. (4^)) should be used. This requirement of high accuracy is 
a consequence of the strong hierarchy among yi and exactly vanishing ee-entry of m. 

For the solutions fi (the first column in eq. ( |4.2[ )) the dominant ^ur-block of m is 
reflected in a dominant 23-block of /. The solutions fi (the second column in eq. ( |4.2| )) 
exhibit a hierarchical structure, with a dominant 33-entry; the maximal 2-3 mixing angle 
in m is generated from a small mixing angle in fi . All the solutions fi and fi have one order 
one eigenvalue, so that for the heaviest RH neutrino mass one finds M3 ~ vr. The other 
two eigenvalues can be as small as ~ 10^^. Type II (I) seesaw contributions dominate the 
entries mg^, and m^r for the solutions /i^2 (/i,2)- The other entries receive significant 
contributions from both seesaw types (hybrid seesaw). This is also true for all the entries 
of the dual pairs of solutions /s, /s and f^, f^. 

2. The largest uncertainty in the structure of nii, is due to the unknown ab- 
solute mass scale of the light neutrinos. Therefore, it is interesting to study the de- 
pendence of the structure of / on this scale. Let us consider the same set of input 
parameters as above, but change the eigenvalues of rrii, as follows: wl(— 0.1, 0.2, 1) — > 
(—mi, ^ ml + Am^^p \/mf + Am^^^). The solutions for / will now depend on mi, which 
may vary between zero and ~ 0.23 eV (due to the cosmological upper bound m^ < 0.7 
eV). The structures of / given in eq. ( |4.2| ) correspond to mi = O.Ivl 0.005 eV. In fig. |5|, 
we plot the masses Mi^2,3 of the RH neutrinos, given by the eigenvalues of Mr = vnf ■, as 
functions of mi. We chose for illustration the first and the third pairs of dual solutions. 

For the solution /i, the mass spectrum of RH neutrinos spans at most one order of 
magnitude, all three masses being rather close to vr. The two lightest RH neutrinos 
undergo a level crossing at mi ~ 0.035 eV. Such crossing points were already identified in 
[23] for the case of pure type I seesaw. In these points the lepton asymmetry generated 
in the decays of the RH neutrinos can be resonantly enhanced, which may be crucial for 
reproducing the observed value of the baryon asymmetry of the universe [23]. For mi ^0.1 
eV (quasi-degenerate light neutrinos), the three RH neutrino masses also become quasi- 
degenerate, since fi becomes a solution with dominant type II seesaw, so that /i ~ m. 

The dual solution fi has a substantially different spectrum. The RH neutrino masses 
span about three orders of magnitude, and two of them can be much smaller than vr. 
Level crossings occur between the two heaviest and between the two lightest RH neutrinos. 
For mi ^0.1 eV, by duality /i corresponds to the dominant type I seesaw. The masses 
Mi of RH neutrinos then scale approximately quadratically with the Dirac-type Yukawa 
couplings Ui. 

For the solution /s, two RH neutrino masses are close to vr for all values of mi, while 
the lightest RH neutrino mass Mi decreases by almost three orders of magnitude when mi 

*Here and below in this section we adopt the numbering of solutions according to which fi go to the 
corresponding solutions given in eq. (4.2), when the input parameters approach those chosen for eq. (4.2). 
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Figure 5: The masses Mi_2,3 of the three RH neutrinos versus the hght neutrino mass scale mi. 
In the upper (lower) panel solid curves correspond to the solution /i (/s) and dashed curves to 
the dual solution /i (/s). We chose the tri-bi-maximal mixing, the light neutrino mass spectrum 
(-mi, ^m^ + Am^^;, y/ml + AmltJ, vl ~ 0.051 eV, vr w 5.9 • lO" GeV and yi.2,3 = lO'^.-i-O. 
The irregularities around mi = 0.011 eV in the lower panel are a numerical artifact. 



increases in its allowed range. For the dual solution, the RH neutrino masses lie between 
10^^ and 10^^ GeV. This dual pair corresponds to hybrid seesaw for all values of mi. 

3. The ordering (normal or inverted) is another important information about the light 
neutrino mass spectrum that is presently missing. Let us consider the input parameters 
chosen in example 1, but replace the normally ordered spectrum of m with an inverted one: 
(-0.1,0.2,1) ^ 
replace eq. (OF) by 



1 + e, 1 + e, e), where e = (l/4)Am^^;/Am^j^ < 0.01. This amounts to 



m 



I 



-2 4 -4 
4 1 -1 
-4 -1 1 



+ e 



/l o\ 

1 
yOOly 



(4.3) 
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Figure 6: The masses Mi_2,3 of the three RH neutrinos versus the hght neutrino mass scale 1713. 
We chose the same input parameters as in fig. but the inverted light neutrino mass spectrum 



\/ml + Am^t„j + m^). In the left (right) panel, the sohd 



curves correspond to the solution /i (/a) and the dashed curves to its dual solution /i (/s). 



Then, one pair of dual solutions is given by 

/-0.33 0.67 -0.67 \ 
... 0.19-0.08 
1.08 



fi 



V 



fl 



I 



( 4 • 10-5 

V ••• 



-0.001 0.001 \ 
-0.01 -0.08 
... -0.91 J 



(4.4) 



and the pair /2, /2 exhibits a qualitatively similar structure. Type II (I) seesaw dominates 
the entries of the first row of f\ (/i), while hybrid seesaw determines the entries of the 
23-block. A third pair of dual solutions is 



V 



-0.35 0.51 -0.52 
... -0.77 0.68 
-1.67 



v 



0.02 0.15 -0.15 
... 0.95 -0.84 
1.85 



(4.5) 



and the pair /4, /4 exhibits a qualitatively similar structure. All entries of these matrices 
receive significant contribution from both type I and type II seesaw. 

In fig. |6| the RH neutrino spectrum is shown as a function of the absolute mass scale 
7713 of light neutrinos. The eigenvalues of my are chosen as {—\ 'rn'^ + Am; 



atm 



m| -|- Am^^„ -|- Am^^;/2, ma), which generalizes the choice fL(— 1 + 1 + e> e) that leads 
to eqs. (4.3)-(^^). For the solution /i, the masses of the three RH neutrinos are rather 



close to each other (within a factor 5) and to vr. For the solution /i, instead, they are 
spread over 3 — 4 orders of magnitude, with the lightest one being around 10^^ GeV. These 
features of the mass spectrum of RH neutrinos are rather insensitive to variations of 
between zero and 0.23 eV. For both /a and /s and in the whole allowed range of m^, the 
masses of two RH neutrinos are close to vr, while the third mass is significantly smaller. 
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Figure 7: The masses Mi. 2. 3 of the three RH neutrinos as functions of sin^ia. We chose the same 
input parameters as in fig. but fixed mi — 0.005 eV and allowed non-zero 613. In the left (right) 
panel, the solid curves correspond to the solution /i (/a) and the dashed curves to its dual solution 

A ifs). 

Actually, despite the fact that all the entries of /a and /a in eq. (|]^) are of order one, the 
determinants A3 and A3 turn out to be much smaller than one, corresponding to one RH 
neutrino mass being much smaller than Vf(. Notice also that the value of M3 is almost the 
same for these two dual solutions. Contrary to the case of the normal mass ordering, the 
present set of input parameters does not lead to the level crossing phenomenon. 

4. Up to now we were assuming in this section the exact tri-bi-maximal mixing. How- 
ever, the present uncertainties of the values of the leptonic mixing angles are quite sizable. 
In particular, only an upper bound on the 1-3 mixing exists, sin0i3 < 0.2. Measuring this 
parameter is one of the main goals of the near-future experimental neutrino program. Let 
us discuss the modifications to the structures of / when a non- vanishing ^13 is allowed. 

Consider the same input parameters as in example 1, but with sin0i3 = 0.1. Then 
eq. ( |4.1D is replaced by 

^ 0.01 0.17 -0.03^ 

0.17 0.53 0.44 . (4.6) 
-0.03 0.44 0.56 j 

The solutions for / are modified slightly with respect to those in eq. ([4.2|), but their 
qualitative features remain the same. However, the relative size of type I and II seesaw 
contributions to a given element of / may change significantly with ^13. The mass spectrum 
of RH neutrinos in general depends weakly on ^13, as shown in fig. |^ for the pairs of dual 
solutions /i^3 and /i^a: the masses Mj change at most by about a factor of 10 for sin^ia 
varying between and 0.2. This indicates that similar underlying theories at the seesaw 
scale may result in very different values of ^13. 

5. An important experimental and theoretical issue is whether CP is violated in the 
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Figure 8: The masses Mi^2,3 of the three RH neutrinos versus the Majorana-type CP violating 
phase p. In the left (right) panel, the sohd curves correspond to the solution /i (/2) and the dashed 
curves to its dual solution /i (/2). We chose the same input parameters as in fig. ^ but assumed 
the quasi-degenerate mass spectrum of light neutrinos ^^[e^^*'', 1 + /S.m'l^i/ {2v\), 1 + Am^j„j/ (2i;|^)] 
with VL = 0.2 eV. 



leptonic sector. The Dirac-type CP-violating phase is associated to the parameter ^13 and 
therefore its effect on the structure of rrii, is small. On the contrary, the two Majorana-type 
CP-violating phases, that is, the relative complex phases of the eigenvalues 772-1,2,3 of f^v^ 
can have large effects on the structure of the mass matrix. In particular, the relative phase 
p between mi and m2 affects substantially the parameter niee = |(7?i.jy)ee|, which determines 
the decay rate of nuclei undergoing neutrinoless 2/?-decay. Since next generation 2/30j/ decay 
experiments will be mostly sensitive to the quasi-degenerate neutrino mass spectrum, we 
consider here the dependence on p assuming three light neutrinos with mass ~ 0.2 eV.^ 

For definiteness, we once again use the parameters chosen in example 1, except that 
the mass spectrum of light neutrinos is now taken to be Ui(e~^*'', 1 -|- e, 1 -|- 77), with 
vl = 0.2 eV, e ~ /S.m'j^^J {2v\) and r] ~ Am^j^/(2?j|^). The mass matrix of light neutrinos 
takes the form 



m 



-2ip 



6 



-4 2 -2 
2-11 
-2 1 -1 



+ 



1 



-1 



1 

1 -1 
-1 1 



+ 



/o o\ 

1 1 
yoily 



(4.7) 



In particular, one finds mee ~ vl\/^ — (8/9) sin^ p. For p = {p = it/2) the CP parity of 

z^i is equal (opposite) to that of the other two light neutrinos. 

The eight solutions for / depend on the value of the phase p in the interval [0, vr] . 

The dependence of the masses of RH neutrinos on p for two pairs of dual solutions is 

^In the basis where y is real and diagonal, the matrix m in the seesaw formula depends on 6 CP-violating 
phases, the three low energy ones plus three phases that are not observable at low energies, but may be 
relevant for the reconstruction of /. They can be taken into account by allowing yi to be complex. In the 
examples considered here we set them equal to zero. 
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Figure 9: The masses Mi^2,3 of the three RH neutrinos as functions of 2/2- In the left (right) panel 
the solid curves correspond to the solution /i (/a) and the dashed curves to its dual solution /i 
if 3)- We chose the same input parameters as in fig. ^, but fixed mi — 0.005 eV and varied 2/2 
between yi = 0.01 and 2/3 = 1. 



presented in fig. ^. The figure is invariant under the transformation p ^ ir — p because, 
in the absence of the other complex phases, this transformation is equivalent to taking the 
complex conjugate of m and therefore of /, so that the masses of RH neutrinos remain 
unchanged. The dependence of Mi on p is generically weak (they change by less than one 
order of magnitude). The three RH neutrino masses corresponding to the solution fi are 
quasi-degenerate for p = (type H seesaw dominates: /i m ~ I3), whereas for p = Tr/2 
they are split by an order of magnitude. For the solution /i (dominant type I seesaw) the 
mass spectrum of RH neutrinos is strongly hierarchical independently of the value of p. 
For the solution /2, level crossing between M2 and M3 occurs at p = 7r/2. 

6. Next, we discuss the dependence of / on the values of the Dirac-type neutrino 
Yukawa couplings yi. In the previous examples we assumed the hierarchical values ys = 
1; 2/2 = 0.1, yi = 0.01. However, as pointed out in section the flavor structure of 
the Dirac-type Yukawa couplings of neutrinos may be qualitatively different from that for 
charged fermions. 

Let us keep fixed 2/3 yi and vary y2 in between. In fig. ^ we plot Mi as functions of 
2/2 for the dual pairs /i^3 and /i^3, taking for the other input parameters the same values as 
in example 1. The left side of the figure corresponds to the limit 2/2 = 2/1 = 0.01, the right 
side to the limit 2/2 = ys = 1- The central value 2/2 = 0.1 corresponds to example 1 (more 



precisely, to the dual pairs /i^s and /i^3 in eq. (|4.2D , that is, to the value mi = 0.005 eV in 
fig. ^). One can see that for the solutions /i^3 the variations of Mj with 2/2 are small; this is 
because type I contribution proportional to 2/2 is subdominant. In contrast to this, for the 
solutions /i^3 the two lightest RH neutrino masses decrease significantly with decreasing 
2/2- For the solution fi a level crossing occurs at 2/2 ~ 0.3. 
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Let us focus on the special case 1/2 = Us, which may be motivated by the experimental 



observation {v^ + z^t-)/\/2, that is, ^13 
the structures of /i and /i : 



and 023 ~ 7r/4. For illustration, we show 



v 



-0.0003 0.14 -0.14 
... 0.86 0.76 
0.86 



/ 0.0003 



-0.04 0.04 
-0.31 -0.31 
... -0.31 



(4.8) 



In this as well as in the other three dual pairs, one finds / 
is a consequence of the choice m, 



12 = -/i3 and /22 = /33- This 
-TUeT, mfj.f_i = rrirT and 7/2 = Us, as can be directly 
verified using eq. ( |3.25|) . Moreover, this "2-3 symmetry" also connects the dual pairs of 
solutions to each other. This can be best seen by considering the seesaw formula after the 
maximal 2-3 rotation: 



nip 



V2m, 



m 



MM 



( 



m 








I fn V2f: 



\ 



12 

/22 — /23 

/22 + /23 J 



(4.9) 



yi 
. . . 2/2 
y2 J 



\ 



yi 

... y2 
... y2 J 



^ hi V2/12 \ Y 

• • • /22 — /23 
\ j \ /22 + /23 y \ 

The equation for the 33-entry decouples from the rest of the system and is quadratic in 
/22 + /23- One its root is the third eigenvalue common to four solutions for /. The other 
root is the third eigenvalue common to the four duals. The equation for the 1-2 block has 
the structure of the LR symmetric seesaw in the case of two generations. As a consequence, 
there are four solutions for /n, /12 and /22 — /23- Each of them can be combined with one 
or the other solution for /22 + /33) giving two solutions which have the same first and second 
eigenvalues. Note that, despite the presence of this "2-3 symmetry" in the neutrino sector 
in the considered example, charged leptons break this symmetry badly because 7^ m,-. 

7. Finally, let us examine the dependence of the masses of the RH neutrinos on 
the dominance parameter d defined in eq. ( p. IT] ). For the one- and two-generation cases 
we illustrated this dependence in figs. |2| and ^, respectively. We now fix all the input 
parameters as in example 1, but allow x to be different from 1. This gives d ~ \?>/x ~ 8-10^^ 
GeV /vr. Notice that eqs. ( |3.25|) and (|3.29| ) imply that the dependence of / on x amounts 
to a common rescaling of the three Dirac-type Yukawa couplings: f{x,yi) = f{l,yi/^/x). 

In fig. IC we plot the masses Mj of the three RH neutrinos as functions of d for the 
dual pairs of solutions /i^3 and /i^a. Since Mj are given by the eigenvalues of vjif, they 
generally decrease as vji (x 1/d with increasing d. However, this is not the case when type 
I seesaw dominates, since in this case / oc 1/vr and Mj tend to constant values, as it was 
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Figure 10: The masses Mi 2.3 of the three RH neutrinos as functions of the dominance parameter 
\d\ = \Av'^y1/{m^)'^^\vL/vji. We chose the same input parameters as in fig. ||, but fixed mi = 0.005 
eV and allowed to vary, so that d — 7.8 • 10^^ GeY /vr. In the left (right) panel, the solid curves 
correspond to the solution /i (/a) and the dashed curves to its dual /i (/a). 

already shown in [13]. Such type I seesaw dominance is realized for d <C 1 for the solution 
/i (left panel in fig. p!o| ). Correspondingly, type II seesaw dominates in /i, so that Mj grow 
linearly with vr and the mass spectrum of the RH neutrinos is proportional to that of light 
neutrinos. As discussed at the end of section |3^ (see also Appendix the other three 
dual pairs exhibit intermediate features for d <C 1, with only Mi {Mi and M2) tending 
to a constant for the solutions /2,3,4 (/2,3,4)) as shown in the right panel of fig. Large 
d corresponds, instead, to a cancellation between type I and II seesaw contributions to 
rriy, that is, ^ —fi, so that the mass spectra of the RH neutrinos for each pair of dual 



solutions tend to coincide, as can be seen on the right sides of the plots in fig. IC. Further 
discussion of the dependence of the mass spectrum of RH neutrinos on vr can be found in 
[13]. 

4.1 Right-handed neutrinos at TeV scale 

The seesaw mechanism is not accessible to direct experimental tests unless the mass scale 
of the seesaw particles is as low as ~ TeV. Therefore, from the phenomenological point 
of view, it is interesting to investigate what regions in the space of the input parameters 
could lead to new particles with masses at this scale. In this section we present a number 
of examples where such a scenario is realized. 

(1) One vanishing or very small Dirac-type neutrino Yukawa coupling (yi 0). This 



limit was studied analytically in section |3.3.l| , where we showed that in this case four 
of the solutions for / have one eigenvalue vanishing as y^. We illustrate this phenomenon 
numerically in fig. 0, where 2/3 and 1/2 have fixed values (1 and 0.01, respectively), whereas 
yi is varied. As expected, in each pair of dual solutions there is one with Mi decreasing 
as yf. The other RH neutrino masses are practically independent of yi when it becomes 
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Figure 11: The masses Mi_2,3 of the three RH neutrinos as functions of yi for two dual pairs of 
solutions. We chose 2/2 = 0.01, nii = 0.005 eV, sin^?i3 = 0.1, while the other input parameters are 
the same as in fig. |^. In both panels dashed curves correspond to the solution with Mi vanishing 
in the limit j/i — > and solid curves to the dual solution. 



smaller than ~ 10~^. Notice that for yi ~ 10~^ the mass of the lightest RH neutrino A'^i 
is already as small as ~ 10 TeV. However, its Yukawa couplings to the SM leptons are also 
suppressed as yi and therefore there is no hope to detect A''i through these couplings. For 
yet smaller values of 7/1, A'^i could be in principle so light as to play a role in cosmology or 
in active-sterile neutrino oscillations. 

(2) Lowering of the LR symmetry breaking scale vr down to ~ TeV. This is allowed 



by direct searches of RH currents, as discussed in section 2.1. Such a scenario is not viable 
for generic ("natural") values of the Yukawa couplings, since it would result in the masses 
of light neutrinos which are much larger than the observed ones. Nonetheless, in view of 
the unique phenomenological possibility to directly test the breaking of the LR symmetry 
and the seesaw mechanism, it is worth accepting the necessary tuning of Yukawa couplings 
and investigate its consequences. 

The simplest phenomenologically viable scenario with vr ~TeV corresponds to choos- 
ing tiny values for all the Dirac-type neutrino Yukawa couplings yi. This allows one to 
avoid conflict with the low-energy neutrino data while maintaining the perturbative uni- 
tarity of the solutions for /, fij < 1. Notice that the seesaw formula is invariant with 
respect to a simultaneous decrease of yi and x provided that yf/x remains unchanged. 
Consider the same choice of the input parameters as in example 1 of section except 
that the values x = 1, 2/3 = 1 are replaced with x = 10"^*^, y^ = 10~^, keeping and 
the hierarchy between yi unchanged. Then the solutions for / are unchanged, but since 
v/j ~ 6 • lO'* GeV in this case, the masses of RH neutrinos are smaller than those in example 
1 of section Q by a factor lO^''. As one can see in fig. ^, this means that the lightest RH 
neutrino can have a mass ^ 100 GeV. Such light RH neutrinos would not be observable 
through their Dirac-type Yukawa couplings, which are tiny; however, because of the LR 
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mixing, they may interact strongly with the ordinary matter through the Yukawa couphngs 
/ with scalar triplets as well as through their coupling to RH gauge bosons. More generally, 
the RH gauge bosons as well as the non-standard Higgs particles related to the LR sym- 
metry breaking would provide clear phenomenological signatures if light enough. If level 
crossing occurs between two TeV-scale RH neutrinos, they could still lead to a successful 
leptogenesis via the resonant enhancement of the produced lepton asymmetry (see section 

13). 

(3) It should be noted that very small values of x, as in scenario 2 above, would indicate 
an approximate symmetry of the scalar potential, since for a generic choice of order one 
couplings one expects vlvr ~ f^, that is x ~ 1. Very interestingly, it is possible to realize 
a TeV scale LR symmetry breaking even preserving this "naturalness" relation. Indeed, 
consider a modification of scenario 2 in which is increased by a large factor W and, 
at the same time, all yf are reduced by W. Then the solutions / of the seesaw formula 
maintain the same matrix structure, but are also reduced by an overall factor W. This 
means, in particular, that the RH neutrino masses are much smaller than vji. Using the 
values of the parameters above and taking W = 10^, one finds Mi ~ keV. Summarizing, 
if LR symmetry breaking occurs at scales much lower than the Grand Unification scale 
and the natural relation vlvji ~ is preserved, this would indicate that the Dirac-type 
and triplet Yukawa couplings y and / are all much smaller than one, with the related 
consequences for phenomenology. 

(4) Finally, let us investigate if TeV scale RH neutrinos may be compatible with at 
least some order one Dirac-type Yukawa couplings and vr ~ TeV. This would require a 
strong fine-tuning in the structure of the Yukawa coupling matrices, in order to cancel 
large contributions to the masses of light neutrinos. To provide a simple example, it is 
convenient to start in the basis where the light neutrino mass matrix is diagonal, nidiag = 
ip^ujdiag/vL = diag{mi, m2, TTis) . Let us define = mi/m2 and choose 

^r^s irs ^ 
y= irs - s . (4.10) 
V ysj 

In this case the seesaw equation for / has only two solutions, the reason being as follows. 
The matrix y in eq. ( [4. 101 ) has one zero eigenvalue, which reduces the number of solutions 
for / from eight to four (section p. 3.1 ). Furthermore, in the basis where y is diagonal, one 



has ruee = 0, which makes two of the four remaining solutions singular and thus unphysical 
(see eq. ( p.34| )). The leftover two solutions are 



xm^i 



(4.11) 
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Notice that the entries of y proportional to s cancel exactly in the seesaw formula and 
therefore do not contribute to the light neutrino masses. As a consequence, the coupling s 
can be as large as one. In contrast to this, y^/x has to be smaller than one to guarantee 
/ss ^ 1- To ensure perturbative unitarity of (/i,2)m, vl should be < eV; therefore a low 
LR scale vr ~ TeV requires < 10~^. 

To consider the structure of the seesaw mechanism in this scenario in the flavor basis, 
it is sufficient to rotate on the left and on the right by the leptonic mixing matrix U: 



m 



U*mdiagU^ , y ^ U*yU^ , /i,2 ^ C/7i,2f/^ ■ (4.12) 



In this flavor basis, type I contributions to the elements of the light neutrino mass matrix 



ruafs are proportional to U*^U^^y'^/x. 



The masses of RH neutrinos are given by the eigenvalues of /i_2 (i.e. the diagonal entries 
in eq. ( 4.12| )) multiplied by vr ~ TeV. The two possible solutions for / are distinguished by 



the mass of the third RH neutrino, while Mi^2 are uniquely determined, with M| — = 
Aiti'^^iVr/vl- Two of the RH neutrinos (A^i,2) have order one Yukawa couplings to the 
light neutrinos and charged leptons. 



5. Structure of seesaw in specific left-riglit symmetric models 

In this section, we discuss the basic features of models which incorporate type I+II seesaw 
mechanism. In particular, we analyze their Yukawa sector and derive the seesaw formula 



for the light neutrino mass in the form given in eq. (1.3). A crucial issue is the structure 
of the matrix of Dirac-type Yukawa couplings of neutrinos y, which is an input in the 
bottom-up approach we adopted to reconstruct /. 

5.1 Minimal LR symmetric model 

We begin by reviewing the structure of the minimal LR model based on the gauge group 
SU{3)c X SU{2)l X SU{2)r x U{1)b-l- It contains the following (color singlet) Higgs 
multiplets: <^(2,2,0), Ai(3, 1,— 2) and A/j(l,3,2), where in the brackets we indicate the 
SU{2)l X SU{2)r X U{1)b-l quantum numbers. The leptons are assigned to L(2, 1, —1) = 
{t^L II)^ and L'^(l, 2, 1) = {N'f^ 11)'^ , and their Yukawa couplings can be written as 

-Cy = ^L'^iaiCALL + ^ L""^ ia2C ArL^ + gL^Ca2^a2L'^ + hL^C^*L^ + h.c. , (5.1) 

Zi Zi 

where a2 is the isospin Pauli matrix and the following conventions are assumed: 

/a+/V^ A++ \ /Ai/V2 A^- \ (4.l4.t\ ,„, 

^^=1 AO -Ai/V2)' A", -A„/V^)' *=(*,-^,J' 
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Denoting the VEVs of the neutral components of ^ as vi^2, ^or the Dirac-type mass matrices 
of the charged leptons and neutrinos one finds 



mi = vue = vig + v^h , uid = vy = v^h + V2g ■ (5.3) 

At low energies the Dirac-type neutrino Yukawa coupling can be rewritten in the form 
of eq. (|L2|) by defining (j)^ = {vi4'f' + Uj^^)/^; where v = -x/luiP + |u2p « 174 GeV is 
the electroweak symmetry breaking parameter. ^'^ Notice that if the theory is supersym- 
metrized, the coupling h should be removed from eq. (^^) due to the requirement of the 
analyticity of the superpotential. As a consequence, y = ye{y2/vi) = ?/etan/3 and the two 
Higgs doublets, (f)u = {4>f 4>2)'^ ^^'^ 4>d = (0? 02 )"^) do not mix.^-*^ 

The Majorana mass matrices of neutrinos are generated by the VEVs fL,_R of the 
neutral components of A^^i^: 

"T-L = VLh , Mr = vr/r . (5.4) 

In the limit where the eigenvalues of of AIr are much larger than those of mo, the mass 
matrix of light neutrinos takes the form 

m,, c^niL- moM'^m^ = vlJl vIr^V^ , (5.5) 

vr 

which is known as type I+II seesaw formula. 

At this point, a discrete LR symmetry should be introduced if one wants to guarantee 
the equality of SU{2)l and SU{2)r gauge couplings, that is, the asymptotic LR symmetry 
of the model. There are essentially two possibilities: the discrete transformation can act 
as a charge conjugation or as a parity transformation. In the previous sections we were 
implicitly assuming the former possibility (we postpone the discussion of the latter one to 
Appendix ly). Then, the discrete LR symmetry acts on matter and Higgs fields as follows: 

L^L", ^>^$^, Al^Ar. (5.6) 

This yields 

f^fL = fR, g = g^, h = h^. (5.7) 

In this case the mass matrix of light neutrinos has the form given in eq. (|1.3| ). Moreover, 



eq. (5/7) implies that y is symmetric, so that the seesaw duality is realized. (If, instead. 



^°The field 0° is thus identified with the neutral component of the SM Higgs boson. The orthogonal 
combination of the neutral scalar fields, 77'^ = (— W2<^i* + Vicj)2)/''^y has zero VEV and so does not contribute 
to the fermion masses. Since it mediates flavour changing neutral currents, it must be heavy and so 
decouples from the low-energy dynamics. 



Similarly to eq. (5^), for quarks one has rud = vyd = viQq + V2hq and = vyu = Vihq + V2gq, where 



the matrices gq and hq are in general independent of g and h, and hq is absent in the SUSY case. 
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the model contained a Higgs bidoublet <!>' with the transformation law <!>' ^ —^'^ un- 
der the discrete symmetry, the corresponding Dirac-type Yukawa coupling y' would be 
antisymmetric, leading to the seesaw duality as well). 

We have seen that in the minimal SUSY LR model the matrix y is proportional to 
the Yukawa coupling matrix of charged leptons ye. Therefore, in the basis where ye is 
diagonal, also y is, and the ratios of its eigenvalues are me/m^ and rrifj^/mr- This scenario 
is ideal for the reconstruction of /, since the input related to the structure of y is completely 
determined. However, one should keep in mind that this minimal model does not provide a 
realistic description of the quark sector, as it leads to cx rud. As follows from eq. 
in the non-SUSY case, the structure of y can be different from ye- For a recent realistic 
model, see [24]. 

Let us briefly consider an upgrade of the minimal LR model to the Pati-Salam gauge 
group SU{2)l X SU{2)r x S'C/(4)c, with a Higgs field in the (2,2, 1) representation. In this 
case quark- lepton unification occurs, which for the Yukawa sector amounts to identifying 
the quark and lepton coupling matrices: Qq = g and hq = h. This implies, in particular, 
y = Uu and ye = Vd- In this case the eigenvalues of y are determined by the up- type quark 
masses, while the mixing is given by the CKM matrix. In non-minimal Pati-Salam models, 
also Higgs fields transforming as (2, 2, 15) may couple to fermions. Their contribution to 
the lepton masses differs from those for the quarks by the factor —3 because of the different 
respective values of {B — L). If both bidoublet types exist and contribute significantly to 
the fermion masses, the direct connection between y and yu is lost. 

5.2 50(10) models 

A natural grand unified embedding of LR symmetric models is provided by 5*0(10) models. 
As is well-known, 50(10) accommodates all the standard model fermions as well as RH 
neutrinos in a single multiplet 16^?. Fermion bilinears transform under 50(10) as 

16f X 16i. = 10, + 120a + 126, , (5.8) 

where the subscript s{a) indicates symmetry (antisymmetry) in the flavor indexes. As a 
consequence, there are three possible types of Higgs multiplets which may contribute to 
fermion masses at the renormalizable level: IOh and 126/^ give symmetric contributions, 
while 120/f gives an antisymmetric one. 

The Higgs triplets ^l,r are both contained in 126//, so that the Majorana mass 
matrices of left-handed and RH neutrinos originate from the coupling 

/ 16f16fT26// 3 f{LLAL + L'L'Ar) . (5.9) 

One generator of 50(10), known as D-parity, acts as in eq. ( |5.6| ). In other words, the 
discrete LR symmetry discussed above is promoted to an automatic gauge symmetry and 
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this guarantees the proportionaUty of the two terms on the r.h.s. of eq. (|5.9| ). The Dirac- 
type neutrino Yukawa couphng y may receive contributions from the Higgs multiplets of all 
three types, provided that they acquire VEVs in the appropriate direction. In fact, both 
IOh and 120/^ contain one bidoublet $ ~ (2,2, 1) under the Pati-Salam group, and both 



120h and 126h contain one (2,2, 15) bidoublet. If only 10^ and 126// contribute to y, 
it is symmetric; if only 120// does, it is antisymmetric. 



Consider first the minimal possibility where only one 126// multiplet develops a VEV in 
the i^L^i direction, that is, y oc /. In this case both type I and type II seesaw contributions 
to nil, are proportional to /. The latter is uniquely determined by mu up to an overall 
factor. There is no seesaw duality in that case, since y is not independent of / and so it is 
not invariant under f ^ m — f . 



If two (or more) 126//i multiplets are introduced, the proportionality between m/, = 
VLifi + VL2f2 and M/j = vmfi + VR2f2 is in general lost. The matrix y depends on the 
same Yukawa couplings fi^2'- V = — 3(u"26 i/i + ^126 2/2); where ^"26 i VEV of the 

"up-type" Higgs doublet (i.e. of the doublet with non-zero VEV in the up-type isospin 



direction) contained in the (2,2,15) component of 126//j. Nonetheless, if vpi2 = 0, an 
effective seesaw duality is realized, with /i and /2 playing the role of / and y, respectively: 

rh = h-Ufi^y. (5.10) 

where 

mu - {VL2 - 18t^r26,2<26,l/^^m) /2 , VliVri - (3^5*26,1)^ - our 

^ = TTH ^27 ' ^ = 2 ' = -3^126,2/2 • 

VLl - {3vl26,l) /VRI 

(5.11) 

If all the parameters in eq. ( p.lO| ) except /i were known, one would be able to solve ( |5.1C| ) 
for /i using the techniques developed in the previous sections for solving the usual seesaw 
formula ( |l.3| ) for /. For three lepton generation this would yield four dual pairs of solutions. 



Consider now the scenario with one 10// and one 126// multiplets. Since 10 is a real 
representation, 10^ transform as 10// under 50(10), so that in general two new Yukawa 
couplings should be added to eq. (^.91) : 



5 16f16/.10// + /i 16/^16/710^. (5.12) 

One then finds 

mD = vy = v'^o9 + vf*oh-3v'^2ef, (5-13) 

where v^^ are the VEVs of the up- and down-type Higgs doublets contained in the (2, 2, 1) 
component of 10//. Since y depends on /, one cannot regard it as an input for solving the 
seesaw formula for /. However, one can consider as an input the matrix vy = v^qQ + vf^h; 
note that in this class of models the couplings to 10// usually give the dominant contribution 
to the mass matrices of the charged fermions. 



-34- 



Then, the seesaw formula (|l.3| ) can be written as 
m = f--yf y, m = , x = ^ . (5.14) 

X VL- {^^^26) /VR 

This equation can be solved for / as usual, taking m, y and x as input parameters. In some 
models of this class, the GUT symmetry breaking is such that only the Higgs doublets in 
IOh develop VEVs (a mixing between these doublets and those in 126h is induced only 
in the presence of 210h multiplets). In this = 0, and one recovers m = m, y = y 

and X = X. Moreover, one has y = yu, that is, the matrix y is fully determined if we 
know the Yukawa couplings of up-type quarks. If also ^^26 ~ 0, one finds ye = yd, which 
fails to reproduce correctly the masses of the down-type fermions of the second and first 
generations. The corrections coming from the 126// (and therefore proportional to vf2Qf) 
may cure this problem. In this case, however, the form of / is constrained not only by the 
seesaw formula, but also by the values of the masses of charged fermions. 

A comment on the number of the Yukawa couplings to 10// multiplets is in order. 
If one wants to study minimal models with only one such Yukawa coupling matrix, there 
are two options: (i) if only one 10// is introduced and the model is supersymmetric, the 
coupling h in eq. (|5.12 ) is automatically forbidden by the analyticity of the superpotential; 
(ii) one can forbid h also in the non-supersymmetric case by assuming that the multiplet 
10// is real. Then only the coupling g contributes to the fermion masses, leading to 
y = yu oc ye = yd, which needs substantial corrections in order to match the differences 
among the mass spectra of charged leptons, up- and down-type quarks. This turns out to 



be extremely constraining, if one tries to fit all the data by adding just one 126// multiplet 
(see [6] and references therein) . If one is willing to consider models with one more Yukawa 
coupling matrix, in the SUSY case one needs to introduce a second 10// multiplet (see e.g. 
[13]), whereas in the non-SUSY case one complex 10// is sufficient. In this framework the 
up sector is decoupled from the down sector, but the relation ye = yd persists. 



Alternatively, instead of adding 10// to the 126// multiplet, one can explore the pos- 
sibility where 120// is added. The couplings in ( [5.12 ) are then replaced by 

g 16f16f120h + h 16f16f120*h , (5.15) 

where the matrices g and h are antisymmetric. Denoting f "2o\ the VEVs of the up (down)- 
type Higgs doublets contained in the (2, 2, 1) component of 120//, and similarly w^io^s 
the (2,2, 15) component of 120//, one finds 

mD^vy= «2o,i - 3<20,i5)5 + (^m,i " ^<20,i5)h - 3<26/ ^yy- 3<26/ , (5-16) 
where y is antisymmetric. Analogously to eq. ( 5.14| ), the seesaw relation takes the form 

r , ~ - VLVR - (3t'i26)^ 

m = f + -yf y, m = — , x = ^ . (5.17) 

X VL - (3^126) /^R ^ 
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This equation can be solved for /, taking m, y and x as input parameters. Using the results 
of Appendix which apply to the case of antisymmetric y, one finds two solutions for /. 

Finally, let us discuss models with both IOh and 120// multiplets added to 126// (for 
an early analysis of fermion masses in this case, see [7]). In the search for the minimal 
realistic SUSY S'O(IO) model, it has been recently shown [25] that the mass matrices of 
the charged fermions can be reproduced through 10// and 120// Yukawa couplings only, 
at least in the first approximation (a three generation fit presents some problems [26]). 
In the case when the coupling / to 126// is also introduced, satisfactory fits of charged 
fermion and neutrino parameters have been recently obtained [27]. In this scenario there is 
no seesaw duality (at least in the form considered here), since y contains both symmetric 
and antisymmetric contributions. Nonetheless, the non-linearity of the seesaw formula still 
leads to multiple solutions for /. To take this into account, one should (i) fix the 10// 
and 120// Yukawa couplings by fitting (approximately) the masses of charged fermions; 
(ii) use the seesaw formula to derive the different structures of / that reproduce a given 
set of neutrino data; (iii) test which solutions for / provide the small corrections needed 
to achieve a satisfactory fit of the masses of charged fermions. 



6. Further considerations 

In this section we briefly discuss the following issues pertaining to our analysis: stability 
of our results with respect to the renormalization group evolution effects and baryogenesis 
via leptogenesis. 

6.1 Stability of the seesaw formula 



Up to now we were assuming that the seesaw formula (1.3) describes accurately the mass 
of light neutrinos at low energy scales. However, due to the breaking of the discrete LR 
symmetry at a scale w//?, the renormalization group (RG) evolution effects below this scale 
can result in a violation of the conditions in eq. (p.7|), which in turn would modify the 
seesaw formula ( |1.3D . A full study of the RG effects in the type I+II LR symmetric seesaw 
is beyond the scope of this paper; here we constrain ourselves to a classification of possible 
effects and estimate of their size. 

Above the scale vlr the discrete LR symmetry ensures fi = fn and y = y^ . Below 
this scale, the fields related by the LR symmetry acquire in general different masses, so 
that // and //? may evolve differently and y may get asymmetric corrections. For example, 
if vlr is larger than vr (the scale of S'C/(2)/j breaking), logarithmic in vr/vlr radiative 
corrections arise. The size of this effect depends on the details of the mass spectra of RH 



^^In 50(10) models, the discrete LR symmetry can be actually broken already at Grand Unification 
scale, even if SU{2)r is not [28]. 



-36- 



gauge bosons and Higgs particles responsible for the breaking of SU{2)fi, which are highly 
model-dependent. We will therefore not discuss this case any further and will just assume 



At this point the fate of the seesaw formula resides with the particle spectrum below 
Vfi. If the masses of all three RH neutrinos as well as of A/, are close to vr, they can be 
integrated out all together, leading to the low-energy effective neutrino mass matrix rriu of 



much smaller than vr, they will contribute to the RG evolution of fi, /r and y down to 
their mass scale. Generically, these corrections induce a splitting between and fn and 
an asymmetry in y. One could envisage two possible approaches to determine the effective 
mass matrix rui, at the lightest seesaw scale Mg (the smallest among Mi and Ma^)- One 
option is to evolve /l,_r and y from vr down to Mg and then integrate out the RH neutrinos 
and Al all together. This has the advantage of determining /l,_r and y as they enter in the 
seesaw formula ( |5.5D . Another option is to integrate out Njn and each at its own mass 
scale, redefining iteratively the effective rriu. Given the /3-functions to some finite order in 
perturbation theory, this second approach should provide a more accurate result. 

For type I seesaw models, a detailed analysis of RG effects, including the running 
between the different mass scales of RH neutrinos, has been performed in [29]. For type 
I-l-H seesaw, such an analysis is not yet available. Here we just estimate the order of 
magnitude of possible effects. One-loop corrections to the matrix elements of / and y can 
be schematically written as 



which, for judicious choices of parameters, should be at or below the percent level. For 
comparison, the present precision of most of the input parameters in rriu is at about 10%. 
Moreover, the loops induced by Njn exchanges are typically proportional to y^. The lighter 
Nfii, the smaller the corresponding Yukawa coupling, since rrip ~ y'^/Mju. Therefore, for 
lighter RH neutrinos, a RG evolution over a wider range may be partly compensated by 
smaller couplings. Similarly, loops induced hy Al exchanges are proportional to and, 
since rrii, ~ v^f ~ v'^ f /M^^, the lighter the A/,, the smaller the expected coupling / and 
the corresponding radiative correction. As a simplified numerical test, we chose certain 
structures of / and y at vr and modified their elements in the seesaw formula by order 
percent corrections. Next, we implemented the usual bottom-up procedure to reconstruct 
/ and found that, among the 8 dual solutions, one reproduces the original structure of / 
within percent errors. However, this is not the case when strong hierarchies among the 
entries of / and y are present and small matrix elements receive corrections proportional 
to the large ones. We expect that only in these special regions of the parameter space 



that the conditions fi = /r and y'^ 



y still hold at the scale vr = vlr- 



eq. (1.3). If, on the contrary, one or more RH neutrinos and/or A^ have masses that are 
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instabilities may occur and large departures of the reconstructed structures of / from the 
true ones may arise. 

To complete the discussion of RG effects, one should consider the evolution of rriu from 
Ms down to the to electroweak scale, where it can be compared with experimental data. If 
only the standard model (or MSSM) particles contribute to the RG evolution, the running 
of the mass squared differences and mixing angles is rather small and negligible for the 
purposes of this paper, except perhaps in the case of quasi-degenerate light neutrinos of 
same CP parity. This could affect the reconstruction of / in example 5 of section 

6.2 Leptogenesis 

As we demonstrated above, in the realistic case of three lepton generations there are eight 
different matrices fi which, for a given y, result in exactly the same mass matrix of light 
neutrinos niu- A natural question is then how one can discriminate between these eight 
possible solutions. The seesaw formula ( |1.3| ) cannot tell us more than it already did, and 
a new independent source of information is necessary. Such a source could be provided by 
the ability of fi to reproduce the observed baryon asymmetry of the universe. 

It is well known that the seesaw mechanism not only explains nicely the smallness of 
neutrino mass, but also has a built-in mechanism for generating the baryon asymmetry 
of the universe through leptogenesis [4]. In this mechanism, first a lepton asymmetry is 
produced in out-of-equilibrium CP-violating decays of heavy RH neutrinos and/or Higgs 
triplets, which is then reprocessed into baryon asymmetry by electroweak sphalerons. The 
amount of the produced lepton asymmetry depends on (i) mass matrix of RH neutrinos 
Mji (masses, mixing angles and CP phases); (ii) Majorana-type Yukawa coupling of leptons 
to Higgs triplets; (iii) Dirac-type Yukawa coupling of leptons to Higgs doublets. In this 
paper, we consider theories with fi = Jr = f, and so the parameters involved in (i) 
and (ii) essentially coincide. This renders the computation of the lepton asymmetry more 
predictive than in generic type I-I-II seesaw scenarios. 

The impact of multiple solutions for / on leptogenesis has been already analyzed in 
a class of 50(10) models in [13] and very non-trivial results were found. Here we will 
not undertake a quantitative analysis of leptogenesis, but rather will make some general 
remarks, which may be of guidance for the development of specific models: 

• It is recognized by now [20, 30, 31] that the presence of both RH neutrinos and 
Higgs triplets may lead to leptogenesis scenarios that are qualitatively different from 
those in pure type I seesaw models. We would like to stress, as a minimal possibility, 
that leptogenesis can work in models with and just one RH neutrino species A'^^ 
(this effective pattern occurs, e.g., when the other two RH neutrinos are super-heavy 
and/or very weakly coupled). Since Nji couples to a unique linear combination L 
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of the flavor eigenstates of lepton doublets, this scenario may be called "flavorless 
leptogenesis" . This is a viable possibility since the relevant interactions of Nji and 
Al contain an unremovable CP-violating phase (which manifests itself in the one- 
generation seesaw formula, see section 3.1). A recent study of this scenario can be 
found in [32]. 

• In many realistic cases an extra freedom gained by the interplay of type I and type 
II seesaw terms may turn out to be insufficient to cure the shortcomings of thermal 
leptogenesis in the pure type I scenario. In particular, it remains true [30, 33] that the 
decaying particle should generically be very heavy (above > 10^ GeV) to produce a 
sufficient asymmetry, thus requiring large reheating temperature with the associated 
problems. Also, the hierarchical structure of Yukawa couplings, which is natural in 
unified scenarios, may strongly suppress the asymmetry. The existence of various 
solutions for /, in particular of those with a non-hierarchical structure, may alleviate 
this problem. Further details can be found in [13]. 

• Given the rich structure of the Higgs sector at the seesaw scale in LR symmetric and 
50(10) models, sources of lepton asymmetry other than the decays of RH neutrinos 
and/or Al into left-handed leptons should not be overlooked. As an example [34], 
a successful leptogenesis may be due to decays of RH neutrinos into a RH charged 
lepton and an SU{2)l singlet charged scalar (contained, e.g., in the 120h multiplet 
of 50(10)). The coupling fji between the RH leptons and Aji (see eq. (|5.lD ) may 
also affect the evolution of the asymmetry. 

• In minimal models with "natural" values of the seesaw scale and of Yukawa couplings, 
in most regions of parameter space the produced lepton asymmetry tends to be too 
small. A way out may be provided by the resonant enhancement of the lepton 
asymmetry which occurs when two RH neutrinos are quasi-degenerate in mass [35]. 
As we demonstrated in several examples presented in section ^, the crossing of the RH 
neutrino mass levels indeed occurs in the LR symmetric seesaw for certain choices of 
input parameters and for some of the solutions . The requirement of level crossing 
may be a powerful constraint for model building. 

7. Discussion and summary 

In models with more than one source of neutrino mass, to disentangle different contributions 
using only low-energy data is not in general possible. The situation is different when these 
contributions have a common origin, as it is the case in left-right symmetric seesaw models. 
In a wide class of such models and their partially unified or Grand Unified extensions, the 
mass matrix of light neutrinos rriiy contains type I and type II seesaw contributions which 
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depend on the same Majorana-like triplet Yukawa coupling matrix /. In this paper we 
undertook for the first time a thorough phenomenological bottom-up analysis of such a 
scenario. We have shown that the interplay of type I and type II seesaw terms in rriu 
may change our interpretation of neutrino data and provide an insight into the underlying 
theory at the seesaw scale. 

We have adopted an approach in which the mass matrix of light neutrinos and the 
matrix of Dirac-type Yukawa couplings y are considered known, and the seesaw relation 
is solved for the matrix of the Majorana-like Yukawa couplings /, which coincides (up 
to a constant factor) with the mass matrix of RH neutrinos and deeply characterizes the 
structure of the underlying theory. To this end, we have developed a linearization procedure 
which allowed us, for symmetric or antisymmetric y, to solve the seesaw non-linear matrix 
equation and obtain exact analytic expressions for the matrix / in a compact form. 

For symmetric y, the overall number of solutions was shown to be 2" for n lepton gen- 
erations. Thus, in the realistic case of three generations, there are eight different matrices 
/ which, for a given y, result in exactly the same mass matrix of light neutrinos rrii,. 

We have studied implications of an intriguing duality property [12] of the LR symmetric 
seesaw mechanism with symmetric or antisymmetric y, which relates pairwisc different 
solutions for the matrix /. The eight solutions of the seesaw equation for / form four 
dual pairs. We have demonstrated that, if one of the solutions / of the seesaw equation 
corresponds to type I dominance, then the dual solution corresponds to type II seesaw 
dominance and vice versa. The other solutions then in general correspond to a hybrid 
seesaw with type I and type II contributions to rrii, being of the same order. An important 
consequence of this result is that, knowing only the values of the input parameters, one can 
determine if there are solutions with one seesaw type dominance, but cannot decide which 
particular seesaw type dominates. 

We have explored the behaviour of the analytic solutions for / in a number of interest- 
ing limiting cases: (a) one seesaw type dominates; (b) the two seesaw contributions almost 
cancel each other; (c) the eigenvalues of y are strongly hierarchical; (d) one light neutrino 
mass vanishes. By analyzing several numerical examples, we found that the masses of 
RH neutrinos exhibit the following generic features: (i) differ strongly for different dual 
solutions for /; (ii) depend crucially on the light neutrino mass spectrum; (iii) depend 
weakly on variations of ^13 within its allowed range and on variations of the low-energy 
leptonic CP-violating phases; (iv) depend strongly on the hierarchy among the eigenvalues 
of y only in the cases with predominant type I seesaw contribution to rrii,; (v) may become 
quasi-degenerate in several regions of the parameter space which are close to the points 
where the level crossings occur; (vi) can be at TeV scale or lighter even when is much 
larger and/or some entries of the neutrino Yukawa coupling matrix y are of order one. 
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seesaw formula 


neutrino Yukawa y 


section 


one gen. 


two gen. 


three gen. 


eq. ( 


O) 


generic (no duality) 


3.1 and 3.2 


2 


4 




eq. ( 


O) 


symmetric 


2.3 


2 


4 


8 


eq. ( 


O) 


antisymmetric 


App. A 


1 


2 


2 


eq. ( 


C.3 


) 


generic (no duality) 


App. C 


2 or 






eq. ( 


C.3 


) 


Hermitian 


App. C 


2 


4 


8 


eq. ( 


C.3 


) 


anti-Hermitian 


App. C 


2 or 


4 or 


8 or 



Table 1: Multiplicity of solutions for the Majorana-type Yukawa coupling matrix /. The first 
column indicates the seesaw equation under consideration, the second column specifies the assump- 
tion made on the structure of the Dirac-type Yukawa coupling matrix y, the third column indicates 
the section where the multiplicity was derived. Fourth, fifth and sixth columns give the number of 
solutions in the case of one, two and three lepton generations, respectively. The dash indicates the 
cases for which the number of solutions was not derived. 

We have also studied (in Appendix ^) the case of antisymmetric Dirac-type Yukawa 
coupling matrices y and showed that our linearization procedure works in that case as well, 
leading to simple analytic expressions for /. The multiplicity of solutions in that case was 
found to be 1, 2 and 2 for one, two and three lepton generations. 

In addition, we considered (in Appendix ^ an alternative realization of the discrete 
LR symmetry, in which it acts as parity transformation rather than charge conjugation. 
This leads to type II and type I contributions to rrii, depending on / and /* rather than on 
the same matrix /. In that case exactly the same duality of solutions holds, provided that 
the matrix y is Hermitian or anti-Hermitian and invertible. In the present paper we gave 
explicit formulas for the solutions for one and two lepton generations, and briefly outlined 
the approach to the three-generation case. The multiplicity of solutions for n lepton gen- 
erations was shown to be 2" in the case of Hermitian y, whereas for anti-Hermitian y there 
are either 2" or no solutions. The latter possibility is a consequence of the non-analytic 
dependence of the seesaw relation on / in this realization of the discrete LR symmetry. 
Our results on the numbers of solutions for / in various cases are summarized in Table 1. 

Since one of our main goals was to provide a guidance for building specific models 
incorporating the LR-symmetric seesaw, we have discussed the structure of the seesaw 
formula in the minimal LR model, Pati-Salam model and in several 50(10) models with 
renormalizable Yukawa couplings. Our analysis shows that in most cases the phenomenon 
of seesaw duality is realized: one can identify a matrix of couplings / which enters in the 
seesaw formula both directly and through its inverse. Such matrix can be reconstructed 
by making use of the methods developed in this paper. 

We have discussed briefly the stability of our results with respect to the renormaliza- 
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tion group running effects, which below the LR symmetry breaking scale may result in 
modifications of the relations fi, = /if, and = ity that were at the basis of our anal- 
ysis. Conditions were found under which these renormalization effects will not spoil our 
reconstruction of the matrix / from the input data. A comprehensive study of this issue, 
however, will require a dedicated effort. 

Finally, we pointed out that the usual leptogenesis mechanism of the generation of 
baryon asymmetry of the universe produces an asymmetry that strongly depends on the 
adopted solution for / and therefore may in principle help discriminate between the eight 
allowed solutions. Moreover, the level crossing between the masses of the two lightest RH 
neutrinos, which occurs for a number of solutions, may provide a resonant enhancement of 
the produced asymmetry. 
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A. Reconstruction of / in the case of antisymmetric y 

Consider the case of antisymmetric Dirac-type neutrino Yukawa coupling, y = —y'^. A cru- 
cial difference with respect to the case of symmetric y is the fact that (2n -|- l)-dimensional 
antisymmetric matrices are not invertible, so that the usual duality among the solutions 
for / does not hold for an odd number of lepton generations. In the case of one generation, 
one trivially has y = 0, and the seesaw is purely of type II. 

The case of two generations is more interesting. An antisymmetric 2x2 matrix is 
defined by a single parameter. Consider the system (|3.9D for y^2 = ^tS = and y^s = 



—yr2 ^ y- Using the same linearization procedure as in section 3.2, one easily finds the 
general solution: 

xX 

f = — zTrm with {xXf - (2y'^ + xdetm)xX + y^ = 0. (A.l) 

xX — y^ 

The solutions X± of the quadratic equation correspond to a pair of dual solutions f±. 
The duality relation ( 3.12 ) is replaced by x^A+A_ = y^. Notice that both solutions are 



proportional to the matrix m. 

While in general the system ( |3.9D has four solutions, two of them are singular when y 
is exactly antisymmetric, as assumed above. This can be understood as follows: consider 
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the case where the diagonal entries of y are equal to zero and the off-diagonal ones satisfy 
UfiZ ~ y ^-nd yr2 ~ —y, so that y is nearly antisymmetric, Then / is still (approximately) 
given by eq. ([A.l|), while the quartic equation for A reads 



[(xA)^ - {2y^ + X det m) xA + y^] (xA - w , 



(A.2) 



to be compared with eq. ( [A.l| ). The pair of dual solutions corresponding to xA — > is not 
physical in the limit, since all the matrix elements of / diverge. 

Consider now the relative size of type I and II seesaw contributions to nii^. The two 
solutions of the characteristic equation for A in eq. ( |A.1| ) can be written as 



xA± 



X det m 



l + ^±(l + d)V2' 



d 



X det m 



(A.3) 



A straightforward calculation shows that both m^^ and ml are proportional to and 
their ratio is 



{mi)i 



l + ^±(l + d)V2 



(A.4) 



For |d| <C 1, r^^^^ ~ —{A/d)^^, so that one seesaw type dominates. For \d\ ~ 1, r^J.^^ ~ 1, 



corresponding to hybrid seesaw. For |fi| » 1, r_|_ 



//// 



-1, which implies a considerable 



cancellation of the two seesaw contributions. These three cases resemble closely the corre- 
sponding limits in the one lepton generation scenario (section |3.1D . 

Let us now turn to the realistic case of three lepton generations. As was pointed out 
above, the usual duality does not apply since y has one vanishing eigenvalue. However, a 
simple analytic reconstruction of the matrix / is still possible, as we show below. 

Any antisymmetric 3x3 matrix can be written as 



yz y2 
-y3 yi 
\-y2 -yi J 



U 



\ 





-2/123 

\ yi23 J 



u 



T 



Uy'U^ . 



(A.5) 



where 7/123 = + + |?/3p and the unitary matrix U is given by 



U 



1 



^ ^1^23 



yi23y23 



-yh 







-^2^23 -^2^1 2/32/123 
\ 2/32/23 2/32/1 2/22/123 J 



(A.6) 



with ^23 = \/|2/2p + |2/3p- Defining m' = U^mll* and /' = U'^fU*, one can write the 
seesaw equation as 



X 



(A.7) 



-43 - 



The entries of the first row of /' are uniquely determined as = m'^j, whereas to find the 
2-3 block one needs to implement the usual linearization procedure. This yields 

_ xXm'22 - yi23m'l2 fi _ xXm'23 - 1/123^2^3 _ xAm^g - yl22.m% 

J 22 ~ \ 2 / ' •'23 ~ \ 2 / ' 33 ~ \ 2 / ' 

~ Uua^ll ~ yi23"^ll ~ yi23"^ll 

(A.8) 

where A is the solution of the quadratic equation 

(xXf - {2yl23m[^ + x det m')xX + y^23"^ n = • (A.9) 

The two solutions X± have the same form as in eq. ( [A. 3D , but with d = 4y^23^n/(^ det m'). 
Even though the usual duality property does not hold in this case, a different duality 
between the two solutions is present: one finds xA+A- = yi23"i'ii and /+ + /i = where 
rhij = m'- + m'^^m'^-/m'ii. The two solutions for / in the original basis are obtained as 
f± = Ufj.U'^, where U is given in eq. ( [A.6| ). 

Thus, for an antisymmetric y the realistic case of three lepton generations can be 
reduced to that of an effective two-generation system, so that there are only two solutions 
of the LR seesaw formula for /, which are related by a modified duality. 



Note that the method of counting the solutions developed in section 2.3 does not apply 
to the case of antisymmetric y: for one and three generations, because y is not invertible, 
and for two generations because this is a degenerate case. Indeed, for n = 2 one finds 
f^^ = [— det det /)] •/, so that the equation m = f + f^^, which replaces eq. ( 2.10| ) in 
the case y'^ = —y, is linear rather than quadratic in /. This equation results in a quadratic 
equation for det/, leading to two solutions for /. 



B. Relative size of type I and II seesaw in the presence of flavor mixing 

In the simple case of one or more unmixed generations the relative size of type I and type 



II seesaw contribution to m,y was discussed in section 3T. Here we analyze this issue in 
the case of 2 and 3 generations with flavor mixing. 

Due to the mixing, each entry of rriiy receives contributions from type I and II seesaw in 
different proportions. For simplicity, we will present only the conditions for the dominance 
of one seesaw type in all the matrix elements of m (extensions to more general cases can 
be easily obtained). In what follows, the relations involving 171^/3 and yi will therefore be 
assumed to hold for each a, /? and i, if not otherwise stated. 



Consider first the two-generation case. From eq. ( 3.10 ) one can see that, when a 
solution Ai satisfies \xXi\ 3> \yiyj\ {i,j = 2,3), one obtains /i ~ m (assuming that m^^ 
and rriT-T are of the same order). This corresponds to the dominant type II seesaw. Then 
the dual solution XA2 = xAi = ^22/3/(^-^1) modulus <C \yiyj\ and the corresponding 
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matrix /2 = /i takes the form obtained in type I seesaw case. In analogy with eq. (^.2]), 
let us define the set of parameters 



i,j = 2,3, a,l3,-i,5 = n,T . 



(B.l) 



Notice that eqs. ( |3.13D and (3.14) imply that a necessary condition to have |xAi| ^ linyjl 
is \mai3m^s\ ^ MUiUj/xl, that is, 

Therefore, this limit ensures the existence of a pair of dual solutions /i^2 with the dominance 
of one seesaw type, in analogy with the limit |d| ^ 1 in the one-generation case. The 
solutions Ai^2 are defined by choosing in eq. ( |3.13 ) and can be expanded in terms of 
the small parameters as 



xAi = X det m 



1 + 



k 



xX2 = xXi 



x(det m) 



x det m 



^ 3(detm)^y|y(^ ^ 
x2(det m)^ 

k 

+ ... 



a;(det m)^ 



(B.3) 



The solutions for / are related by f\ 



m 



/2 with 



1 



X det m 



( 



\ 



-vl 



rrir 



+ 



uln^fifi det m — mT-rk 
rc(det m)^ 



1 



2/22/3 det m + /c 



-vl 



rc(det m)^ 



y2mrT det m 



x(det m)^ 



(B.4) 

The second terms in the square brackets represent the leading order correction to the pure 
type I seesaw. 

In general, condition ( |B.2| ) does not guarantee that the other dual pair of solutions 
/3,4 corresponds to one seesaw type dominance. If no special cancellations occur, eq. (|B.2| ) 
implies, in particular, |detm| ^ 4:\yiyj/x\. If this condition is satisfied, then and /4 
are of hybrid type. This novel feature with respect to the one-generation case is a genuine 
effect of flavor mixing. In fact, the solutions corresponding to the choice of r„ in eq. ( |3.13| ) 
are given, up to higher orders in small parameters, by 

k±s ' 



xXs,' 



2 det m 



It 



a; (det m)^ 



+ ... 



s = VA;2 - (2^22/3 det m)2 . (B.5) 



Notice that A3^4 are dual to each other. The corresponding solutions for / (to leading 
order) are 




(s ± A;) =F 2mT-T-y| det m m^r{s ± A; it 2^22/3 det m) 

mrris ± /c) =F 2m^^y| det m 



(B.6) 
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Roughly, {fz,A)ij ~ f^ap and ^3^4 = A3^4 ~ yiUj/x. Therefore, type II and type I contri- 
butions to are of the same order: m^^ ~ VLiriap and ~ v'^Hiyjfkil {vrF) ~ VLtnap; 
this proves that the solutions 73^4 are of hybrid type. 

If a special cancellation in det m occurs, also the solutions /3^4 are dominated by one 
type of seesaw. Indeed, in the limit det m — > eq. ( |3.13| ) becomes 



1,3 



1 

±- 
2 



kx + 42/22/3 + \kx 



x\ 



2,4 



x\ 



1,3 



2 



kx + %22/3 — \kx 



(B.7) 

where k ~ {jn^^y^ + mT-Ty2)^- When eq. (B^) holds, one has \kx\ S> 4y2y3- As a 
consequence, |xAi| = |xA3| » \yiyj\, so that /i^3 are solutions with dominant type II 
seesaw. By duality, /2,4 correspond to the dominant type I seesaw. 

Note that a suppression in detm-j^ is phenomenologically motivated if the neutrino 
spectrum has normal mass hierarchy. In this case the full 3x3 mass matrix nrii, is dominated 
by large entries in the 2x2 ^r-block. This block has to incorporate the maximal mixing 
observed in the atmospheric neutrino oscillations as well as the hierarchy Am^^^ <C Am^^^. 
These two requirements imply that the determinant of the ^r-block is suppressed. 

When one abandons the condition ( |B.2| ), that is, when at least one jc^^^^^l ^ 1, all 
four solutions are of hybrid type. In other words, some sets of input parameters necessarily 
imply a hybrid seesaw scenario. 



In the limit when all \d]^p^s\ ~ '''±1^ ~ =1=2/22/3) and eq. (3.10) implies 

that (at least the diagonal) matrix elements of / grow as l/(|xA| — |?/22/3|)i thus ending up 
to be much larger than ma(j. This means that one enters the cancellation region, where 
type I and II contributions are almost equal in absolute value and opposite in sign. In 
this limit, the matrix elements of / will violate perturbative unitarity. This means that 
for some values of the input parameters, there are no physically acceptable solutions for /. 
In other words, the LR symmetric seesaw mechanism cannot reproduce certain sets of the 
input parameters and thus, in principle, can he ruled out. 

Consider now the three-generation case. In analogy with eq. ( |B.1| ), one finds that a 
necessary condition for the dominance of one seesaw type is 



VL 








VR 






xruapm^i 



<!, i,j = 1,2,3, a,/3,7,(5 



e,/U,T. 
(B.8) 

In this case, the solutions Aj of eq. ( p. 29 ) can be expanded in the small parameters d^^^^, 
which we denote collectively by d. The largest solution turns out to be 

A 3y2s(detm)2 -^2 



Ai = det m 



1 + 



(B. 



(detm)2 (detm)^ 

Noting that ~ d^, it is straightforward to verify with eq. ( |3.25| ) that fi corresponds 

to the dominant type II seesaw. As a consequence, the three RH neutrino masses satisfy 
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-^1,2,3 ^ yiUjMs, where Mg is the seesaw scale defined as Ms = v'^/m^ (normal mass 
ordering) oiv^/m2 (inverted mass ordering). 

Three more solutions are given by A2,3,4 = A2,3,4[l + 0{d)\, where A2,3,4 are the roots 
of the cubic equation 

A^ + + y25A + y2detm = 0. (B.IO) 
det m 

It is easy to see that A2,3,4/Ai ~ d. From eq. ( |3.25| ) one finds 

(A? - y^detm)mi,- + XlAij - Y'^XkSa 

{fk)ii^— ^ TT—^ k = 2,3,4. B.ll 

^■'^ ^ Xl-Y^XkS -2Y^ det m ' ' ' ^ ^ 

One can verify that for these three solutions type I and II contributions to nii, are of the 
same order (hybrid seesaw). Since Xk = det fk are suppressed with respect to Ai, only two 
RH neutrinos are have masses that are larger than the seesaw scale, M2,3 ^ yiUjMs, while 
Ml ~ ViVjUs. 

The dual solutions, A2,3,4 = — y^/A2,3,4, are suppressed by a factor ~ with respect 
to Ai. The corresponding structures for / are 

^^'^'^ -Xl-Y^-X,S-2Y^deim ' ^ = 2' 3, 4. (B.12) 

Note that fk + fk = m, as required by duality. The solutions fk also correspond to type 
I and II contributions to rrijy being of the same order (hybrid seesaw). The masses of the 
RH neutrinos satisfy M3 » y-iyjMs, Mi^2 ~ Uil/jMs. 
Finally, the solution dual to Ai is given by 

y2 



Ai 



det m 



(B.13) 



This is the smallest solution, since Ai/Ai ~ d"^. By duality, it corresponds to the dominant 
type I seesaw, as one can check explicitly expanding eq. ( |3.25 ). All three RH neutrino 
masses are approximately at the seesaw scale, Mi^2,3 ~ yiVjMs. 

As already mentioned, the parameters d^^^^ may not all be small at the same time. 
In that case, for any of the eight solutions /, only some of the elements of the matrix rrii, 
may be dominated by one seesaw type, while one or more other elements receive significant 
contributions from both seesaw types. 

A remark is in order on the special case when detm is suppressed. For example, one 
may have |detm| ^ \Y\ even if eq. ( |B.8| ) holds. Clearly, in the limit detm — > the 
expansions ( |B.9D -( B.13| ) do not apply. We showed in section 3.3.2 that, when detm = 0, 
two pairs of dual solutions for A are equal in absolute value and opposite in sign to the other 
two dual pairs. Moreover, the equation for A has the same form as in the two- generation 
case (with the redefinitions ( |3.37| )). In particular, the relative size of type I and II seesaw 
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contributions to rrij^ may be evaluated as we did above for the two-generation case: when 
eq. (|B.8| ) holds, there are two opposite solutions itAi which lead to the dominant type II 
seesaw, while their duals correspond to the dominant type I seesaw. The other two opposite 
pairs of dual solutions correspond in general to hybrid seesaw, unless a special suppression 
in \fA occurs. 



C. Alternative realization of the discrete left-right symmetry 

Consider a different realization of the discrete LR symmetry in the minimal LR model, in 
which it acts as parity transformation. Eq. ( |5.6| ) is then replaced by 




L^(LT= / ' $^$t, Al^AIj. (C.l) 



This possibility is almost as popular in the literature (see, e.g., [24, 36]) as the one we 
adopted before, even though it does not allow an 50 (10) embedding of the minimal LR 
model. Eqs. (P-lj) and (|5.1|) imply 



/^/l = /^, 9 = 5^ h = h^. (C.2) 
In this case the seesaw relation is 



m.^vLi -—y[rr'y. (C.3) 

where y is defined in eq. ( [5.31) . The difference in complex phases between the seesaw 
formulas (|1.3| ) and ( |C.3| ) can, in principle, lead to a different phenomenology if CP is 
violated in the leptonic sector. 

It can be easily checked by direct substitution into eq. ( |C.3| ) that, if the matrix y 
is Hermitian or anti-Hermitian and invertible, this equation possesses the same duality 
property as eq. ( |1.3[ ): namely if a matrix / solves eq. ( C.3 ), so does its dual / = rriy/vL — f ■ 



The above conditions on y are both necessary and sufficient for the duality to hold. 

Eqs. ( C.2[ ) and (5.3) imply that y is actually Hermitian if /i = (which holds, e.g., in 



the SUSY case) or if the two VEVs of the bidoublet vi and V2 are both real (this is not 
the case in general, since both VEVs can be non-zero, and their relative phase cannot be 
rotated away [36]). If y is Hermitian, one can diagonalize it according to y = UydU'^ with 
U unitary, and go to the basis where y is diagonal and real by absorbing U in the definition 
of rriy and / (analogously to eq. (|2.8D). 

Consider first the one-generation case with an arbitrary complex y. Defining / = \ f\e^'^, 
one can obtain from eq. (C.3) 

=^<^ - "^1^1 (C.4) 



- y'^/x 
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The value of |/| is determined by the requirement that that numerator and denominator 
on the right hand side have the same moduh (recall that in our conventions all the VEVs 
are real and positive, and so is x): 



i/r 



I ,2 , 2cosx y 
m H 



21 



0. 



(C.5) 



Here we used 
only if 



y\ e^^. Eq. ( C.5 ) has (two) real and positive solutions for |/| if and 



\m\ > (1 — cosxj 



(C.6) 



This means that, if the input parameters do not obey this inequality, no value of / satis- 
fies the seesaw formula ( p.3|) . In other words, such a set of input parameters cannot be 
explained by the LR symmetric seesaw mechanism. This problem never occurs in the case 
of eq. ( |1.3| ), which is analytic in the elements of / and thus always has solutions. 

The requirement that y be Hermitian or anti-Hermitian reduces in the one generation 
case to the condition that y be real {x = 0) or purely imaginary (x = vr). Let us focus on 



the X = case. Eq. (|C.q ) is then always satisfied and the two solutions of eq. (|C.5|) are 



\f±\ = ^ il + \d\)'/'±l 



(C.7) 



where d was defined in eq. (|3.2|) . The discussion of the dependence of on \d\ in 
different limiting cases is analogous to that in section 3A . Finally, the phases (/>± = arg f± 
are determined uniquely by plugging eq. ( |C.7] ) into eq. (C^). Thus, if y is real, eq. ( C.3 ) 
has two solutions. 

Consider now two lepton generations, limiting ourselves to the case y = y^ and working 
in the basis where y is diagonal and real. The linearization procedure is the same as in 



section 3.2. However, to close the system of equations one has to consider both eq. (C.3) 
and its complex conjugate, so that the linearized system contains 6 equations. Defining 
A = |A|e*'', one arrives at the solution which is analogous to that in eq. ( p. 10 ): 



/ 



x|A| 



vlvl 



x|A|m^^ + yim*^^e'P 



(x|A|)2 

The parameter A is determined, as usual, from the equation A = det/(A), which reads: 



(C.8) 



3;3|A|2(|m^^|2y2 + 2\mf,r\'^y2y3 + Imrrpyi) 



^|A| [detme-*^x2|A|2 + {detmye'Py^l] =0. 



(C.9) 



The imaginary part of this equation implies p = arg A = arg det m+lir with Z = 0, 1 (barring 
the special cases detm = and x|A| = y2yz)- Then |A| satisfies a quartic equation, which 
is completely analogous to that in eq. ( 3.11[) . Therefore, its four solutions are easily found 
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analytically and one can verify that for each / = and / = 1 there is only one pair of dual 
solutions which yield real and positive |A| and therefore are acceptable. As a consequence, 
also in this scenario there are four solutions for the matrix /. 

The three-generation case can be considered quite analogously. The linearization pro- 
cedure yields a closed set of 12 linear equations for the elements of the matrices / and /*, 
and solving the characteristic equation one arrives at eight possible solutions for /. 

It is actually not difficult to show that the number of solutions for n lepton generations 
is always 2", provided that the matrix y is Hermitian and invertible. Indeed, going into the 
basis where y is diagonal an real and performing the same transformations as in section 



2.3| , one arrives at the equation 

rh = f-{h-\ (CIO) 

where rh and / are symmetric. The matrix / can be diagonalized as / = VfdV'^, where 
V is unitary and diagonal and real; it is easy to see that then (/*)^^ is diagonalized 
by the same transformation. Therefore, multiplying eq. ( p.lOj ) by on the left and by 
V* on the right, one diagonalizes its right hand side, and so also the left hand side. Thus, 
y is a unitary matrix that diagonalizes the matrix rh and therefore is determined by the 
known quantities m and y. In the diagonal basis eq. ( |C.10| ) yields 

fdi-rndifdi-'^ = 0, i = l,...n. (C.ll) 
This gives two values of fdi for each i: 



/.. = ^±y^+i. (C.12) 

A full ra-generation solution is obtained by picking one of these two values for each i, which 
yields 2"^ solutions for the matrix / and therefore for /. 

The same method can also be applied to the case of anti-Hermitian y (making use of 
the fact that iy is Hermitian). Eq. ( |C.12D is then replaced by 



/.. = ^±y^-i, (C.13) 

This leads to 2" solutions for /, provided that all rhdi satisfy rh'^- > 4; otherwise the seesaw 
equation ( |C.3| ) has no solutions. 

A procedure similar to that described below eq. ( |C.10D (but making use of a complex 
orthogonal transformation rather than of a unitary one) was first used in [13] in order to 
solve eq. ( |1.3D with symmetric y. In both cases, this approach allows not only to count the 
solutions of the seesaw equation, but also to obtain them explicitly. 
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